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SECTION  I 

INTRODUCTION  AND  SUMMARY 


1.  BACKGROUND 

This  document  is  the  final  report  of  contract  F33615-80-C-1218,  NEW  SENSOR  CON¬ 
CEPTS.  The  objective  of  this  effort  was  to  assess  the  feasibility  of  promising  new 
strapdown  inertial  sensors,  using  innovative  approaches.  The  scope  of  the  effort  included 
theoretical  and  experimental  investigations,  laboratory  experiments,  the  building  and 
testing  of  breadboard  accelerometers  and  recommendations  for  further  work  needed  to 
improve  the  potential  of  the  concept. 

This  program  focused  on  developing  an  accelerometer  which  utilizes  force-sensitive 
vibrating  quartz  beams  as  the  basic  sensing  elements.  This  approach  results  in  a  device 
that  is  virtually  a  solid  state  accelerometer  with  an  inherent  digital  output.  The  sensor 
also  has  good  potential  for  low  cost. 

The  Kearfott  approach  in  the  application  of  quartz  crystals  to  accelerometers  has  been  to 
use  flexure-mode  resonators  in  a  beam  configuration.  The  behavior  of  a  vibrating  beam  in 
tension  is  somewhat  like  a  string  in  tension;  an  increase  in  tension  will  cause  the  resonant 
frequency  to  increase.  The  beam,  however,  has  several  advantages  over  the  string: 

•  requires  no  bias  tension 

•  responds  to  compression 

•  frequency  change  with  tension  is  more  linear  than  the  string 

•  achieves  usable  full-scale  force  induced  frequency  change  of  up  to  +10% 
resulting  in  high  instrument  scale  factor 

•  scale  factor  temperature  sensitivity  less  than  10  ppm/°C 

•  predictable  and  hence  modelable  temperature  sensitivity,  an  advantage  for 
heaterless  operation  over  wide  temperature  range. 

In  an  accelerometer  application,  the  crystal  resonator  and  oscillator  circuit  replaces  the 
torquer  coils,  torquer  magnets,  pickoff  and  capture  electronics  of  a  conventional  force 
rebalance  accelerometer.  This  results  in  the  following  advantages: 

•  the  accelerometer  becomes  virtually  solid-state 

•  the  mechanical  assembly  becomes  very  simple  for  low  cost  and  high  reliability 

•  the  electronics  (oscillator)  becomes  very  simple  for  low  cost,  low  power  and 
high  reliability  with  an  inherently  digital  output  (no  pulse  capture  or  A/D 
conversion) 

•  higher  performance-to-cost  ratio  per  acceleration  sensing  channel. 

The  most  successful  applications  of  the  quartz  vibrating  beam  have  been  the  High 
Accuracy  Vibrating  Beam  Accelerometer  (VBA),  for  strategic  missile  and  gravity  meas¬ 
urement  applications,  and  a  commercially  available  pressure  sensor.  In  each  of  the  above 
applications  the  force-sensitive  crystal  technology,  which  is  about  15  years  old,  is 
successfully  competing  against  the  conventional  electromagnetic  force  rebalance  tech¬ 
nology  (D'Arsonval  galvanometer  mechanism)  which  is  about  100  years  old. 

Because  this  force-sensitive  crystal  concept  has  been  very  successful  in  both  ultraprecise 
as  well  as  commercial  use,  its  application  to  tactical  systems,  with  accuracy  requirements 
between  these  two  extremes,  appears  a  natural  choice. 


2.  PERFORMANCE  GOALS  AND  ACCOMPLISHMENT  SUMMARY 

Below  is  a  goal-by-goal  review  of  what  has  been  accomplished  and  what  this  program 
indicates  can  be  accomplished  on  future  units. 

•  Bias  stability:  100  us  Goal  -  Short  term  bias  stabilities  of  less  than  1  n  g  have 
been  demonstrated.  Long  term  stability  is  currently  limited  by  trending,  but 
we  recommended  several  ways  to  eliminate  this  trending  and  achieve  100  ng 
long  term  stability. 

•  Bias  Temperature  Sensitivity;  lMg/°/C  (0°C  to  60°C)  Goal  -  Unmodeled  bias 
temperature  sensitivity  of  2^g/°0  has  been  demonstrated.  Improvements 
along  with  modeling  have  the  potential  of  achieving  l^g/°C  sensitivity. 

•  Scale  Factor  Linearity:  0.001%  Goal  -  Test  data  and  analysis  indicate  that  a 
method  of  matching  and  trimming  to  achieve  the  goal  appears  feasible. 

•  Scale  Factor  Stability:  100  ppm  Goal  -  Analysis  and  test  data  indicate  that 
scale  factor  stability  achieved  is  on  the  order  of  0.1  ppm. 

•  Scale  Factor  Temperature  Sensitivity;  1  ppm/0  C  Goal  -  Scale  factor 
temperature  sensitivity  on  the  order  of  1  ppm/°C  has  been  achieved,  however, 
sensitivities  of  10  ppm/C  are  expected  to  be  more  typical.  This  sensitivity 
will  be  very  repeatable  and  can  be  modeled  and  compensated  to  1  ppm/°C. 

•  Dynamic  Range:  t  10  g  Goal  -  A  +10  g  range  was  easily  achieved;  we 
recommend  that  future  units  have  a  greater  full  scale  input. 

•  Reaction  Time:  1  min  Goal  -  Accelerometer  reaction  is  virtually  instant  on 

with  an  ultimate  goal  of  no  temperature  control. 

•  Scale  Factor  Symmetry:  2  ftg/g  Goal  -  The  approach  taken  is  not  a  pulsed 
device.  Therefore,  this  requirement  is  not  applicable. 

•  Threshold:  50  qg  Goal  -  Thresholds  of  less  than  1  fig  have  been  demonstrated. 
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SECTION  II 

VBA  BASIC  OPERATION 


1.  BASIC  MECHANIZATION 

The  basic  acceleration  sensing  element  in  the  VBA  is  a  pendulously  supported  proof  mass 
constrained  along  its  axis  of  freedom  by  a  crystal  vibrating  beam  force  transducer.  These 
transducers  change  their  natural  frequency  in  a  highly  systematic  way  with  force;  the 
frequency  increasing  with  tension  and  decreasing  with  compression.  The  transducer  is 
used  as  the  resonating  element  in  an  oscillator  circuit.  The  frequency  of  the  oscillator  is 
thereby  a  sensitive  function  of  input  acceleration. 

2.  VIBRATING  BEAM  FORCE  TRANSDUCER  (RESONATOR) 

The  behavior  of  a  vibrating  beam  in  tension  is  described  in  Appendix  A.  The  analysis 
describes  the  relationship  between  beam  geometry,  material  physical  properties,  boundary 
conditions  and  their  effect  on  the  instrument's  bias  frequency  and  frequency-force 
sensitivity.  This  behavior  is  most  conveniently  expressed  as  a  power  series  used  in 
Equations  (1)  through  (3). 

The  quartz  beam  material  has  a  unique  combination  of  mechanical,  chemical  and  thermal 
properties  needed  for  frequency  standards.  In  addition,  quartz  is  piezoelectric  which 
affords  a  very  convenient  way  to  sustain  the  oscillations. 

The  beam  vibrates  in  flexure  as  a  fixed-fixed  beam,  described  analytically  in  Appendix  A, 
and  illustrated  schematically  in  Figure  1. 

The  crystals  are  driven  piezoelec trically  by  means  of  electrode  pairs  arranged  to  provide 
an  electric  field  across  the  minimum  dimension  of  the  beam  and  induce  the  desired 
flexure  stresses. 

The  resonators  themsleves  are  low  cost  items  of  which  thousands  have  been  manufac¬ 
tured.  The  raw  material,  quartz  crystal,  is  relatively  abundant,  in  both  its  natural  and 
synthetic  form,  and  is  used  extensively  in  the  frequency  control  and  communication 
industries.  Machining  techniques  have  been  developed  to  produce  them  in  large  quantities 
at  low  cost.  Many  of  the  machining,  plating,  and  other  manufacturing  processes  used  are 
similar  to  those  used  to  make  the  quartz  crystals  used  in  wristwatches.  Wristwatch 
crystals  also  vibrate  in  flexure  in  a  tuning  fork  or  free-free  beam  configuration.  In  large 
quantities,  they  cost  watch  manufacturers  less  than  two  dollars  each. 

3.  DUAL  BEAM  VBA  MECHANIZATION 

A  schematic  representation  of  the  dual  beam  VBA  is  presented  in  Figure  1. 

The  two  beams  and  two  proof  masses  used  are  arranged  so  that  an  input  acceleration 
places  one  beam  in  tension  and  the  other  in  compression.  The  VBA  output  is  then  taken  as 
the  difference  frequency  between  the  two  beams. 

This  two-beam  mechanization  is  described  analytically  by  Equations  (1)  through  (3).  Note 
that  these  equations  express  the  frequency  behavior  of  the  individual  beams  as  a  power 
series.  The  K  terms  are  the  product  of  the  proof  mass  and  the  appropriate  K  terms  of 
Table  2-7  of  Appendix  A. 


DUAL  BEAM 


FIGURE  1.  SCHEMATIC  OF  DUAL  BEAM  VBA  APPROACH 


BEAM  1  OUTPUT 
fi =  K01  *  Kna  +  K2ia2 
BEAM  2  OUTPUT 


+  K3!a 


f2  =  K02  '  K12a  *  K228  *  K32a 
DIFFERENCE  OUTPUT 


to  =  K„  -  K 


*(K|1*  K12)a  +  (K21-  K22)!2+  (><31  +  K32)a3 


Nominally  Zero 


-Nominally  Double 

1 — Nominally  Zero  -  Common  Mode  Rejection 
a  =  Input  Acceleration 


Nominally  Double 


(1) 

(2) 

(3) 


The  advantages  of  this  dual  beam  push-pull  mechanization  are: 

•  The  difference  frequency  and  hence  accelerometer  bias  is  low.  The  tolerance 

on  the  match  of  the  individual  beam  K  will  determine  the  bias. 

o 

•  The  scale  factor  is  twice  that  of  the  individual  beam  scale  factor. 

o 

•  The  difference  frequency  a  term  is  nominally  zero  to  reduce  vibration  effects 
and  increase  linearity. 
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•  The  difference  frequency  a  term  is  nominally  doubled,  but  is  very  small. 

•  Differencing  results  in  common  mode  rejection  to  greatly  reduce  thermal 
effects. 

The  common  mode  rejection  feature  greatly  reduces  the  effects  of  bias  error  sources  that 
similarly  affect  the  individual  beams  such  as  thermal  effects.  Examples  of  VBA  data 
illustrate  that  the  difference  frequency  is  at  least  an  order  of  magnitude  more  stable  than 
the  individual  beam  frequencies.  This  is  attributed  to  the  beams  being  in  very  close 
physical  proximity  and  therefore  experiencing  the  same  environment.  This  residual 
temperature  effect  on  bias  is  a  result  of  the  mismatch  between  individual  beam 
temperature  characteristics,  kept  small  by  careful  selection  of  beam  pairs. 

The  effectiveness  of  the  common  mode  rejection  principle  is  shown  in  Figure  2.  It  is 
obvious  from  this  example  that  the  randomness  of  the  difference  frequency,  in  equivalent 
input  acceleration,  is  at  least  an  order  of  magnitude  less  than  that  of  the  individual 
beams.  Similar  to  scale  factor,  the  residual  temperature  effect  on  bias  is  very 
predictable  and  therefore  also  modelable  for  heaterless  operation. 
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SECTION  III 

LOW  COST  VBA  DESIGN  APPROACH 

1.  VBA  TECHNOLOGY  BASELINE 

Prior  to  this  program,  Kearfott  had  directed  its  development  efforts  toward  a  high 
accuracy  VBA  instrument  intended  for  strategic  applications.  As  a  result  of  these  efforts, 
Kearfott  now  has  several  capabilities  applicable  to  the  low-cost  VBA: 

•  Force-sensitive  crystal  resonator  design  techniques  were  developed  which 
make  extensive  use  of  finite  element  analysis. 

•  Quartz  crystal  selection,  fabrication,  processing,  plating  and  testing  tech¬ 
niques  were  developed. 

•  VBA  assembly  facilities  and  procedures  were  developed,  including  high  vacuum 
techniques. 

•  VBA  test  procedures  were  developed,  which  included  automatic  data  acquisi¬ 
tion  and  reduction  using  an  HP21MX  Computer  capable  of  monitoring  12 
frequency  signals  simultaneously. 

2.  LOW-COST  TACTICAL  VBA  APPROACH 

Preliminary  experiments  and  investigations  indicated  that  the  dual  beam  VBA  had  good 
potential  for  low-cost  tactical  applications.  To  realize  this  potential,  technology 
development  efforts  were  expended  which: 

•  Increase  sensitivity  by  using  more  of  the  beam  force  measuring  range. 

•  Introduce  a  proof  mass  damping  mechanism  to  protect  the  beams  from  shock 
and  vibration  and  simplify  the  design. 

•  Develop  a  saturating  oscillator  which  has  the  potential  for  low  cost  and  high 
reliability. 

Details  on  development  approaches  that  realized  these  goals  are  described  next. 

3.  INCREASED  SENSITIVITY 

The  VBA,  designed  for  high  accuracy  strategic  applications,  had  no  internal  means  of 
damping  the  beam-proof  mass  mechanism.  Vibration  protection  depends  on  keeping  the 
natural  frequency  of  the  beam-mass  system  high  (>2000  Hz),  and  then  providing  an 
external  vibration  isolator  with  a  lower  (=900  Hz)  break  frequency  to  assure  adequate 
input  vibration  attenuation  at  the  natural  frequency.  Keeping  the  natural  frequency  high 
was  accomplished  by  keeping  the  mass  of  the  proof  mass  low  and  keeping  the  axial 
stiffness  of  the  beam  high. 

Both  of  these  actions  result  in  lower  acceleration  sensitivity  since  a  low  mass  proof  mass 
results  in  lower  forces  and  a  stiff  beam  results  in  lower  frequency-force  sensitivity. 

The  sensitivity  of  the  low-cost  VBA  was  increased  by  increasing  the  mass  of  proof  mass. 
The  actual  amount  of  sensitivity  increase  depends  on  size,  linearity  and  other  considers 
tions,  but  for  the  10  g  full  scale  VBA  developed  under  this  program,  the  factor  was  about 
3.  This  resulted  in  a  16  g  proof  mass. 
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This  increase  in  sensitivity  will  proportionately  reduce  thermal  crystal  aging  and  other 
effects  on  bias.  However,  unless  measures  were  taken,  increasing  sensitivity  also  will 
decrease  the  natural  frequency  and  increase  susceptibility  to  shock  and  vibration.  To 
compensate  for  this,  a  beam  protection  mechanism  and  proof  mass  damping  mechanism 
was  devised. 

4.  PROOF  MASS  DAMPING  MECHANISM 

Proof  mass  damping  was  accomplished  by  a  squeeze-film  gas  gap  between  the  proof  mass 
and  housing.  This  is  a  proven  damping  mechanism  that  has  been  successfully  used  in  many 
accelerometers.  Sufficient  damping  can  be  introduced  to  achieve  a  critical  or  near- 
critical  damping  ratio  and  the  scale  factor  can  be  greatly  increased.  The  resonant 
frequency  of  the  beam-mass  system  can  be  reduced  to  about  300  to  400  Hz.  Enough 
compliance  can  be  built  into  the  resonator  so  that  mechanical  stops  can  be  used  to  protect 
the  beam  from  vibration  and  shock  overload.  This  gas  damping  approach  is  the  most 
simple  and  results  in  lower  cost. 

The  equations  describing  the  action  of  a  gas  damping  mechanism  are  included  as  Appendix 
B.  They  are  for  an  ideal  arrangement  and  do  not  account  for  the  effects  of  nonsymmetry, 
the  central  resonator  mounting  hole  or  the  effects  of  c^  #  c2,  all  shown  schematically  in 
Figure  2;  was  made  longer  than  c^  because  there  is  more  motion  at  the  c„  point.  The 
plan  was  to  assemble  a  unit  and  then  vibrate  it  with  various  gases  and  gas  pressures  to 
empirically  establish  a  performance  baseline.  The  initial  values  selected  were: 

Radius  R  =  1.14  cm 

Gaps  c  =  0.00254  cm,  c„  =  0.00381  cm 
Beam  Spring  Rate  K0  =  5.77  (10 ')  dyne/cm 
Proof  Mass  m  =  16  g 

The  performance  goals  were: 

Natural  frequency  (undamped)  fN  =  300  Hz 
Damping  ratio  =  on  the  order  of  one 

Vibration  test  results  for  various  gases  and  gas  pressures  are  presented  in  Section  4. 

5.  COMPLIANT  RESONATOR 

The  previous  resonator  designs  for  the  high  accuracy  VBA  are  axially  very  stiff  because  of 
the  high  natural  frequency  requirements.  For  the  low-cost  VBA  under  development  an 
entirely  new  resonator  was  needed,  one  that  was  axially  compliant.  The  resonator  design 
arrived  at  is  illustrated  in  Figure  3.  The  only  characteristics  retained  from  the  previous 
resonator  designs  were  the  vibrating  beam  geometry,  the  isolator  structure  and  mounting 
structure,  were  completely  changed.  Because  of  mounting  requirements  and  vibration 
isolation  requirements,  the  resonator  is  a  rather  complicated  structure  with  many  degrees 
of  freedom.  If  not  properly  designed,  the  vibration  frequency  of  one  of  these  degrees  of 
freedom  will  be  close  to  the  fundamental  frequency  of  the  vibrating  beam.  This  can  cause 
activity  dips  which  are  a  well-known  effect  in  the  crystal  industry.  Activity  dips  are 
characterized  by  a  very  erratic  vibrational  frequency  and  a  loss  in  Q.  The  vibration 
isolator  must  also  be  carefully  designed.  The  vibration  isolator  isolates  the  bending 
moment  and  shear  reactions  at  the  root  of  the  beam,  and  prevents  vibrational  energy  from 
escaping  to  the  mounting  surface.  Excessive  escaping  vibrational  energy  will  result  in  a 
decreased  resonator  Q.  In  the  design  of  the  resonator,  extensive  use  is  made  of  finite 


7 


OSCILLATOR 


element  analysis  where  a  candidate  beam  design  is  broken  up  into  more  than  100  separate 
elements.  The  model  is  then  analytically  subjected  to  an  axial  loading  and  the  dimensions 
of  the  isolator  springs  are  adjusted  to  obtain  the  correct  axial  compliance.  Using  this 
procedure,  the  value  of  5.77x10'  dyne/cm  was  arrived  at  in  order  to  obtain  a  300  Hz 
natural  frequency  with  the  16  g  proof  mass.  The  finite  element  model  is  then  subjected  to 
a  mode  frequency  analysis  where  the  vibrational  characteristics  of  the  various  degrees  of 
freedom  are  determined.  An  example  of  the  results  obtained  from  this  analysis  is 
presented  in  Appendix  C.  The  illustrations  presented  in  Appendix  C  are  computer 
generated  graphics,  and  illustrate  the  various  modes  and  frequency  of  vibration  that  cun 
occur  in  the  resonator  structure.  The  mode  and  frequency  predictions  of  this  analysis 
method  have  been  found  to  be  very  helpful  and  very  accurate  in  determining  the  resonator 
vibrational  characteristics.  Observe  that  mode  12  is  that  of  the  fundamental  beam 
vibration  which  is  predicted  to  occur  at  40  kHz  where  indeed  it  actually  does.  The  lower 
frequency  modes  allow  a  check  on  the  behavior  of  the  vibration  isolation  system.  The 
analysis  is  also  used  as  a  criterion  to  be  sure  that  the  various  modes  are  at  least  8  to  10 
kHz  away  from  the  40  kHz  fundamental  to  avoid  an  activity  dip  problem. 

6.  ACCELEROMETER  DESIGN  APPROACH 

Details  on  each  of  the  piece  parts  illustrated  in  Figure  3  follows. 

The  base  serves  as  both  the  bottom  damping  plate  and  the  attachment  point  for  the 
housing  end  of  the  resonator  structure.  It  includes  the  two  headers  which  communicate 
with  the  resonator  electrodes.  Also  included  is  a  copper  exhaust  tube.  Through  this  tube, 
the  device  is  evacuated  during  bake  out  and  back  filled  with  a  predetermined  pressure  of 
gas.  The  design  of  this  piece,  as  well  as  the  other  pieces,  is  very  simple  and  should 
present  no  future  manufacturing  problems. 

The  proof  mass  is  also  illustrated,  although  from  the  view  it  was  photographed  from  the 
flexure  joints  are  not  readily  visible.  Since  the  proof  mass  of  this  VBA,  as  well  as  other 
VBA's,  is  relatively  heavy  compared  to  conventional  accelerometers,  tl.e  flexures  are 
easier  to  fabricate.  This  is  because  a  lapping  procedure  is  not  needed  to  bring  them  down 
to  very  small  dimensions.  The  thinnest  flexure  dimension  used  for  the  proof  mass  was 
0.0012  inches.  The  flexures  were  used  in  the  as-ground  condition  and  did  not  require  any 
additional  lapping.  The  top  and  bottom  surfaces  of  the  proof  mass  serve  as  the  damping- 
surfaces.  The  upper  end  of  the  resonator  is  mounted  the  central  opening  of  this  proof 
mass. 

The  damping  plate  serves  just  that  function. 

As  a  pre-assembly,  the  proof  mass  is  sandwiched  between  the  base  and  the  damping  plate. 
The  gaps  are  set  by  the  circular  shims  shown.  This  three  piece  assembly  is  then  inserted 
into  the  housing,  after  which  a  seal  is  made  using  a  high  temperature  epoxy.  The  units 
assembled  under  this  program  are  essentially  breadboard  units.  To  facilitate  these 
assemblies,  epoxy  was  liberally  used.  The  crystal  was  epoxy  mounted  to  the  assembly,  and 
then  the  base  sealed  to  the  housing  by  a  bead  of  epoxy.  This  approach  resulted  in  easy 
assembly  to  allow  quick  evaluation  of  the  compliant  resonator  and  the  damping  mechan¬ 
ism,  since  these  two  items  were  a  radical  departure  from  previous  high  accuracy  VBA's. 


The  saturating  oscillator  is  also  illustrated.  This  particular  one  was  available  in  flat  pack 
form  and  was  used  during  most  of  the  tests.  It  is  cemented  to  the  bottom  portion  of  the 
base  after  the  exhaust  tube  is  pinched  off  and  tucked  into  its  recess. 


An  example  of  a  complete  assembly  is  illustrated  in  Figure  4.  To  make  a  complete  dual 
beam  VBA  package,  two  of  these  assemblies  are  needed.  The  base  of  each  assembly  has 
recessed  screw  holes  so  that  they  can  be  easily  mounted  back  to  back.  It  is  in  this  back- 
to-back  condition  that  the  devices  were  tested. 

Liberal  use  of  epoxy,  however,  is  a  shortcoming  of  this  present  design  because: 

•  The  ultimate  in  crystal  stability  is  not  obtained  until  a  bake  out  at  least  200°C 
can  be  obtained.  The  epoxy  llimited  this  to  100°C. 

•  The  crystals  being  epoxied  into  place  is  also  a  serious  shortcoming.  Because  of 
temperature  changes  and  changes  in  load  the  epoxy  has  a  somewhat  plastic 
behavior.  This  results  in  creeping  changes  in  the  boundary  conditions  of  the 
beam  and  results  in  the  trends  observed  during  tests. 

For  future  development  units  it  will  be  a  goal  to  have  an  all-welded  assembly.  Means  of 
achieving  this  all-welded  assembly  capable  of  a  full  200°C  bake  out  are  presented  in 
Section  6. 


SECTION  IV 
TEST  RESULTS 


1.  NOMINAL  BIAS  AND  SCALE  FACTOR 

The  nominal  scale  factors  of  the  individual  resonators  of  units  LC1  and  LC3  are  88  and 
89  Hz/g  respectively.  This  results  in  a  difference  frequency  scale  factor  of  nominally 
177  Hz/g. 

The  nominal  bias  frequencies  of  the  individual  resonators  of  units  LC1  and  LC3  are 
40,856.7  and  40,507.2  Hz  respectively.  This  results  in  a  nominal  difference  frequency  of 
349.5  Hz,  which  is  equivalent  to  (after  dividing  by  177)  1.97  g  bias.  It  is  characteristic  of 
the  VBA  to  have  a  high  but  very  stable  bias.  Bias  on  production  units  can  be  reduced  by 
better  matching  of  individual  resonators. 

2.  BIAS  STABILITY 

The  major  deviations  in  design  approach  of  the  low-cost  VBA  over  the  previous  high 
accuracy  VBA  were: 

•  The  introduction  of  gas  damping,  which  means  that  the  resonators,  which 
formerly  operated  in  a  high  vacuum,  would  now  have  to  operate  in  a  fraction 
of  an  atmosphere  of  the  gas  to  be  chosen. 

•  The  former  high  accuracy  unit  used  rather  expensive  linear  oscillators.  The 
low  cost  units  operated  with  relatively  inexpensive  and  simple  saturating  type 
oscillators. 

One  of  the  first  tests  conducted  was  to  compare  the  behavior  of  the  unit  to  these  two  new 
operating  conditions.  The  test  results  are  presented  in  Figure  5.  Figure  5A  is  an  example 
of  short  term  stability  obtained  on  unit  LC4  which  was  operating  under  an  evacuated 
condition  with  a  resonator  Q  of  about  40,000  and  using  one  of  the  relatively  expensive 
linear  oscillators.  Figure  5B  is  a  similar  unit  LC3  operating  in  1/20  of  an  atmosphere  of 
neon  which  resulted  in  a  resonator  Q  of  12,000.  The  oscillator  used  here  was  of  the 
saturating  type.  Very  little,  if  any,  short  term  performance  degradation  was  realized  by 
going  to  these  new  conditions.  Note  that  the  test  results  presented  in  Figure  5  are  for 
individual  accelerometer  halves  tested  separately.  Therefore,  there  is  no  benefit  of 
common  mode  rejection  of  thermal  errors. 

Figure  6  is  an  example  of  bias  stability  obtained  over  a  longer  period  of  time.  Figure  6A 
is  a  plot  showing  how  the  drifts  of  the  individual  resonators  varied  during  the 
approximately  20  hour  period.  As  mentioned  in  Section  II,  the  actual  VBA  output  is  taken 
as  the  difference  frequency  between  these  two  individual  frequencies.  This  difference 
frequency  is  plotted  in  Figure  6B.  Note  that  bias  stabilities  of  well  under  1/ig  are  being 
realized.  Also  note  that  this  is  raw  data  and  has  not  been  compensated  for  any  effects  at 
all.  These  tests  were  run  in  the  +1  g  position  so  that  the  drift  error  effects  illustrated  are 
those  of  both  bias  and  scale  factor.  The  temperature  of  the  room  monitored  during  the 
period  of  the  test  appears  in  Figure  6C.  Because  of  thermal  insulation  and  temperature 
control,  the  resonators  experience  only  about  one-tenth  of  these  room  temperature 
variations.  Note  that  the  periodic  movement  of  the  beam  frequencies  of  Figure  6  is  well 
correlated  with  the  changes  in  room  temperature  of  Figure  6C.  Note,  however,  that  these 
variations  are  absent  from  the  difference  frequency  plot  of  Figure  6B.  This  is  a  typical 
example  of  the  common  mode  rejection  features  obtained  from  the  dual  beam 
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FIGURE  5.  LOW  COST  VBA 

SHORT  TERM  BIAS  STABILITY  INDIVIDUAL  RESONATOR 


design.  The  combination  of  doubling  the  scale  factor  and  common  mode  rejection  has 
improved  the  bias  stability  by  about  one  order  of  magnitude. 

Figure  7  is  a  typical  example  of  a  longer  drift  run.  However,  note  that  the  difference 
frequency  is  trend  corrected  for  this  example.  It  has  been  characteristic  of  the  low-cost 
VBA's  under  test  to  exhibit  first  order  trends  which  are  very  predictable.  The  difference 
frequency  trend  during  this  test  was  about  10^g/day.  It  is  these  trends  that  will  be 
corrected  on  the  future  designs  by  eliminating  the  epoxies  and  changing  to  an  all-welded 
design  so  that  the  crystals  can  receive  benefit  of  200°C  bake  out. 

3.  BIAS  THERMAL  SENSITIVITY 

A  goal  of  the  low-cost  VBA  program  is  to  eventually  provide  an  accelerometer  which  can 
operate  over  a  wide  temperature  range  without  heater  control.  With  respect  to  bias,  the 
dual  beam  approach  depends  on  a  good  thermal  match  between  the  temperature  vs. 
frequency  variation  of  the  two  individual  resonators.  (The  bias  temperature  sensitivity 
model  is  presented  in  Appendix  E.)  Since  the  output  is  taken  as  the  difference  frequency, 
a  perfect  match  will  result  in  a  zero  bias  temperature  sensitivity.  It  has  been  found 
during  the  course  of  oscillator  development  that  a  change  in  the  value  of  the  trim  resistor 
in  the  oscillator  circuit  will  change  the  frequency  of  the  resonator  slightly.  This  is  an 
action  similar  to  "crystal  pulling"  used  to  final  trim  frequency  standard  crystals.  It 
became  apparent  that  if  the  trim  resistor  could  be  replaced  with  a  thermistor-resistor 
network  designed  so  that  the  temperature  sensitivity  of  this  circuit  would  compensate  for 
the  temperature  sensitivity  of  the  resonator,  it  may  be  possible  to  force  a  match  between 
the  two  resonators.  This  experiment  was  tried  and  the  results  are  presented  in  Figure  8. 
Figures  8A,  B  and  C  show  the  before  compensation  results.  This  was  a  dynamic  test 
whereby  the  temperature  change  was  induced  by  simply  turning  off  the  heaters.  As  a 
result,  a  dynamic  thermal  condition  exists.  Figure  8A  and  B  indicate  that  there  is  about 
an  8%  mismatch  in  the  linear  temperature  coefficient  of  the  individual  resonators. 

Figures  8D,  E  and  F  show  the  results  of  the  same  experiment  run  after  the 
thermistor-resistor  temperature  compensation  network  was  substituted  for  the  oscillator 
trim  resistor.  Figure  8E  indicates  that  now  the  temperature  vs.  time  traces  of  the  two 
beams  are  virtually  superimposed. 

A  comparison  of  Figures  8C  and  F  illustrates  that  the  best  straight  line  slope  of 
temperaturg  vs.  bias  has  been  reduced  from  an  original  30jug/°F  down  to  about  lng/°F 
over  the  30°F  temperature  range  of  the  experiment. 

The  results  of  this  experiment  demonstrate  that  very  low  bias  temperature  sensitivity  is 
achievable.  For  production  units,  the  procedure  will  be  even  more  effective  since  a  large 
number  of  low-cost  VBA  "halves"  will  be  available  to  initially  obtain  a  lower  thermal 
mismatch  than  the  8%  of  LC1  and  LC3  breadboard  units. 

4.  SCALE  FACTOR  STABILITY 

It  is  estimated  that  the  scale  factor  stability  is  better  than  1  ppm,  exclusive  of  thermal 
sensitivity  (to  be  discussed  in  Section  IV,  para.  5).  The  following  analysis  is  offered  as  a 
rationale  for  this  estimate. 

To  describe  the  bias  frequency  and  the  force -frequency  effects  of  an  individual  resonator, 
Equation  (4)  is  taken  from  Appendix  A.  Note  that  just  the  bias  and  scale  factor  effects 
are  considered  and  the  terms  are  as  defined  in  Appendix  A. 
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In  Equation  (5),  the  tension  term  is  replaced  by  a  mass  times  acceleration  term  (mg)  and 
Equation  (4)  is  then  rearranged. 


Frequency  Scale  Factor 


Equation  (5)  describes  for  the  individual  resonator  how  both  bias  and  scale  factor  are 
dependent  on  various  characteristics  such  as  beam  boundary  conditions,  beam  geometry 
and  material  properties.  Note  that  both  bias  and  scale  factor  are  dependent  on  these 
various  characteristics  to  approximately  the  same  order  of  magnitude.  That  is,  a  1  ppm 
change  of  any  characteristic  will  change  the  bias  or  the  scale  factor  by  0.5  ppm  to  2  ppm. 
Now,  for  example,  assume  the  aQ  term,  which  is  dependent  on  the  boundary  conditions, 
changes  by  1  ppm  due  to  creep  or  the  epoxy  joint  which  fastens  the  beam  to  either  the 
proof  mass  or  the  base.  This  1  ppm  change  of  a„  results  in  a  1  ppm  change  in  the  bias 
frequency  and  also  a  1  ppm  change  in  scale  factor.  For  the  individual  resonator 
experiencing  1  g,  the  nominal  value  of  the  bias  term  is  40,000  Hz,  while  the  nominal  value 
of  the  scale  factor  term  is  88  Hz.  The  1  ppm  change  of  a«  therefore  causes  a  0.04  Hz 
change  in  the  bias  term  and  just  a  0.000088  Hz  change  in  the  scale  factor  term. 

The  conclusion  indicated  is  that  instability  sources  on  individual  resonators  are  common  to 
both  bias  and  scale  factor  to  about  the  same  degree.  But,  since  the  bias  is  much  greater 
than  the  scale  factor  (about  450:1  at  1  g).  Virtually  all  instabilities  experienced  during 
testing  can  be  attributed  to  the  bias. 
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Because  of  the  performance  limitation  caused  by  epoxy  creep  and  the  inability  to  achieve 
a  high  temperature  bake  out  as  described  in  Section  7,  individual  resonator  frequency 
instabilities  of  about  0.1  ppm  are  being  realized.  This  indicates  that  one  or  a  combination 
of  several  terms  of  Equation  (5)  are  on  the  order  of  0.1  ppm  unstable.  By  the  reasoning 
presented  in  the  previous  paragraph,  this  infers  that  the  scale  factor  is  about  0.1  ppm 
unstable. 

5.  SCALE  FACTOR  TEMPERATURE  SENSITIVITY 

A  series  of  four  position  tests  were  performed  over  the  temperature  range  of  85°F  it 
135°F.  The  dual  beam  VBA  model  equations  for  this  procedure  appear  in  Appendix  D. 
From  experience  on  previous  VBA  designs,  we  normally  expect  a  slope  of  scale  factor  vs. 
temperature  on  the  order  of  five  to  ten  ppm/°F  for  an  individual  resonator.  Because  of 
the  two  beam  mechanization  of  frequency  differencing,  the  VBA  scale  factor  is  nominally 
twice  that  of  the  individual  resonator  scale  factor  and  does  not  receive  the  benefit  of 
common  mode  rejection.  Therefore,  even  though  the  scale  factor  is  double  that  of  the 
individual  resonator,  the  temperature  sensitivity  on  a  ppm/°F  basis  remains  the  same  and 
therefore  we  would  normally  expect  five  to  ten  ppm/°F  scale  factor  temperature 
sensitivity  which  will  be  very  repeatable  and  verv  modelable.  For  this  series  of  tests, 
however,  one  of  the  units,  LC1,  had  a  -7.0  ppm/°F  scale  factor  slope  while  LC2  had  a 
+6.6  ppm/F  slope.  The  reason  why  one  scale  factor  had  a  negative  slope  and  the  other  a 
positive  slope  is  not  currently  known.  The  implication  is  that  for  this  particular  case, 
there  was  a  common  mode  rejection  mechanism  over  this  temperature  range  of  scale 
factor.  Therefore  the  scale  factor  sensitivity  of  the  difference  frequency  was  less  than 
l  ppm/°F.  This  residual  scale  factor  temperature  sensitivity  will  also  be  very  repeatable 
and  very  modelable.  In  production  units,  this  scale  factor  common  mode  rejection  may 
not  be  achieved.  In  any  event,  the  scale  factor  temperature  sensitivity  is  expected  to  be 
no  greater  than  10  ppm/°F. 

6.  NONLINEARITY 

The  same  four  position  test  procedure  of  Appendix  D  also  gives  us  some  indication  of 
nonlinearity,  although  the  £est  conditions  were  rather  crude  due  to  the  bias  instability. 
For  LC1  and  LC3  the  Hz/gZ  coefficients  were  -0.146  and  -0_.^J8  respectively.  For  a  16  g 
proof  mass,  this  translates  to  a  force  sensitivity  of  -5.94  (10  1U)  Hz/dynez  (for  LC1)  and  - 
7.24^0  1U)  Hz^dyneZ  (for  LC3).  Both  these  values  are  somewhat  higher  than  the  -4.44 
(10  10)  Hz/dynez  analytically  predicted  in  Table  2-7  of  Appendix  A.  It  is  not  presently 
known  why  there  was  a  deviation  between  individual  resonators  and  what  has  been 
theoretically  predicted.  However,  the  test  procedure  used  had  shortcomings  because  of 
the  bias  instability,  and  determining  this  coefficient  involves  small  differences  of  large 
numbers. 

However,  it  is  encouraging  that  the  correct  sign  and  order  of  magnitude  was  obtained.  In 
production  units  there  will  be  a  large  number  of  individual  VBA  halves  to  choose  from. 
Also  it  is  anticipated  that  future  VBA's  will  not  have  as  heavy  a  proof  mass.  In  the  dual 
beam  mechanization  the  nonlinearity  terms  do  receive  benefit  of  the  cancellation  so  that 
for  an  8  g  proof  mass  a  VBA  g  coefficient  of  about  4/cg/gz  will  be  achieved  for  a  1% 
mismatch  of  individual  resonator  g  coefficients.  By  trimming,  this  can  be  brought  down 
to  on  the  order  of  lng/g£  or  less.  This  trimming  is  accomplished  very  simply  by  changing 
slightly  the  mass  of  one  of  the  proof  masses.  In  summary  then,  production  units  will 
achieve  high  linearity  by  first  making  a  careful  match  of  the  individual  g  coefficients, 
and  then  trimming  slightly  the  mass  of  one  of  the  proof  masses  to  bring  the  gL  coefficient 
nominally  to  zero. 
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TheQg  coefficients  were  determined  from  multi-position  test  data  over  the  range  of  85  to 
135°F  and  to  the  accuracy  of  the  test  there  was  no  detectable  temperature  sensitivity  of 
the  g  square  coefficient.  This  is  not  surprising  since  the  resonator  theory  presented  in 
Appendix  A  shows  that  the  nonlinear  terms  as  the  scale  factor  and  bias  are  dependent  on 
the  same  characteristics,  such  as  beam  boundary  conchtions,  beam  geometry  and  quartz 
material  properties.  Therefore,  the  sensitivity  of  the  g  coefficient  should  also  be  on  the 
order  of  one  to  ten  ppm/°F. 

7.  ALIGNMENT  STABILITY 

As  a  result  of  the  multi-position  tests  comes  the  alignment  term.  Although  there  was  no 
attempt  to  align  the  devices,  some  insight  was  gained  into  the  type  of  stability  that  we 
could  expect  over  the  temperature  range.  By  nature  of  the  VBA  construction  the 
alignment  is  simply  a  matter  of  the  differential  expansion  between  the  quartz  and  the 
material  of  the  VBAQbody  itself.  Therefore,  an  alignment  sensitivity  of  temperature  on 
the  order  of  1  //rad/  F  will  result  for  each  ppm/°F  mismatch  between  quartz  and  the 
material  body.  It  is  planned  to  make  a  match  of  better  than  1/2  ppm/°F  in  these 
expansion  characteristics.  Therefore  an  alignment  stability  with  temperature  on  the 
order  of  1/2  //rad/F  should  be  expected  on  each  unit.  The  actual  alignment  acheived 
during  the  four  position  test  for  LC1  and  LC3  were  3.8  (10~3)  radian  and  -1.2  (10-4) 
radian,  respectively.  The  test  data  was  not  accurate  enough  to  actually  measure  a  slope 
and  confirm  the  1/2  to  1  /xrad/°F  sensitivity  expected. 

8.  VIBRATION  TESTS 

One  of  the  first  units  assembled  was  subjected  to  vibration  tests,  before  it  was  sealed,  so 
that  some  evaluation  of  the  damping  mechanism  could  be  gained.  The  test  was  run  on  a 
single  unit  so  therefore  there  is  no  benefit  of  common  mode  rejection.  The  object  of  the 
test  was  to  determine  what  the  natural  frequency  and  maximum  transmissibility  was  and 
also  to  determine  if  there  were  any  unusual  resonances.  Figure  9  is  a  good  example  of  the 
type  of  response  that  was  obtained.  The  test  was  run  by  placing  the  output  of  the 
oscillator  into  a  frequency  discriminator  and  then  plotting  the  output  of  the  discriminator. 
Figure  9  shows  that  the  system  behaves  as  a  classical  second  order  system  having  a 
resonant  rise  at  the  natural  frequency.  It  was  a  design  goal  to  place  this  frequency  at  300 
Hz.  It  came  in  at  about  250  Hz.  On  future  designs,  this  frequency  is  easily  altered  by 
adjusting  the  stiffness  of  the  compliant  resonator  or  the  mass  of  the  proof  mass  or  both. 
The  resonant  frequency  was  found  to  be  relatively  insensitive  to  pressure.  This  is 
contrary  to  the  analysis  of  Appendix  B,  which  accounts  for  the  gas  spring  effects. 
Transmissibility  was  greater  (damping  lower)  than  expected  however,  which  agrees  with 
Appendix  B. 

A  series  of  vibration  tests  was  performed  at  various  absolute  pressure  levels  with  various 
gases.  The  results  of  these  tests  are  plotted  in  Figure  10.  The  damping  achieved  is 
relatively  insensitive  to  pressure  down  to  1  to  2  in  Hg  absolute.  As  the  table  in  Figure  10 
shows,  the  value  of  damping  is  proportional  to  the  viscosity  of  the  gas  used.  It  was 
fortunate  to  find  that  very  little  absolute  pressure  is  needed  to  achieve  the  desired 
damping.  This  is  fortunate  because  the  lower  that  gas  pressure,  the  better  the  Q  of  the 
resonators.  This  is  pointed  out  in  Section  IV,  para.  9  where  the  effects  of  gas  on  resonator 
performance  is  discussed.  Because  of  these  test  results,  the  remaining  units  were  all 
filled  with  1.5  in  Hg  absolute  of  neon.  This  limited  the  maximum  transmissibility  to  about 
1.5. 
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FIGURE  9.  LOW  COST  VBA  VIBRATION  TEST 
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Future  development  studies  and  tests  are  recommended  to  determine  the  ideal  gas  or  gas 
mixture  which  will  result  in  both  high  resonator  Q  and  the  desired  damping. 

9.  RESONATOR-GAS  EFFECTS 


Figure  11  summarizes  the  gas  effects  on  resonator  operation.  Figure  11A  is  a  plot  of  gas 
density  vs  viscosity.  As  pointed  out  earlier,  a  gas  of  high  viscosity  provides  the  most 
effective  damping  while  a  gas  of  low  density  results  in  the  least  gas  effects  on  resonator 
operation.  From  Figure  11 A  it  can  be  seen  that  the  gases  which  offer  the  most  favorable 
viscosity  to  density  ratios  are  neon  and  helium.  Figure  11B  is  a  plot  of  resonator  Q  vs  gas 
absolute  pressure.  From  this  plot  two  things  are  apparent.  First,  it  can  be  seen  that 
helium  has  the  least  effect  on  resonator  Q.  Second,  the  relationship  of  Q  to  pressure 
becomes  very  nonlinear  at  low  pressure.  The  reason  for  this  is  not  known,  however,  it  is 
suspected  that  at  low  pressures,  the  mean  free  path  of  the  gas  becomes  commensurate 
with  the  resonator  vibration  amplitude.  This  is  the  extent  to  which  this  phenomena  has 
been  investigated.  Figure  11C  is  a  plot  of  resonator  frequency  change  vs  gas  absolute 
pressure.  This  effect  has  been  found  to  be  linear  and  different  for  the  three  gases 
investigated.  The  pressure-frequency  sensitivity  appears  to  be  approximately  proportional 
to  the  gas  density.  This  observation  is  reasonable  since  it  is  concluded  that  this 
frequency-pressure  sensitivity  is  a  mass  effect.  The  gas  immediately  surrounding  the 
vibrating  beam  is  also  set  into  motion  and  therefore  the  added  mass  of  the  gas  increases 
the  effective  density  of  the  quartz.  It  is  not  surprising,  then,  that  a  gas  such  as  helium 
with  a  very  low  density  would  have  the  least  effect.  The  data  presented  in  Figure  11, 
along  with  the  vibration  test  results  presented  in  Section  IV,  para.  8  indicate  that  the  1/20 
of  an  atmosphere  of  neon  used  may  not  be  the  ideal  condition.  Helium,  or  perhaps  a 
mixture  of  helium  and  neon,  may  be  the  best  choice  to  bring  about  the  best  combination 
of  resonator  performance  (as  indicated  by  high  Q)  and  proof  mass  damping,  (shown  by  low 
transmissibility).  Neon  was  used  in  the  breadboard  accelerometers  because  they  used 
epoxy  seals  and  it  was  feared  that  helium  might  leak  through  these  seals.  However, 
future  designs  will  be  completely  welded  assemblies  which  should  solve  any  possible 
helium  containment  problem.  We  plan  to  do  further  development  work  in  this  area. 
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FIGURE  11.  RESONATOR  GAS  EFFECTS 
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SECTION  V 

OSCILLATOR  ELECTRONICS 


1.  BACKGROUND 

We  have  addressed  the  need  for  a  low  cost,  highly  reliable  oscillator  design  to  complement 
the  dual  beam  low  cost  VBA  approach.  Included  in  this  effort  is  a  review  of  future 
packaging  and  system  interface  requirements  that  will  optimize  cost  and  performance. 
We  will  take  full  advantage  of  the  commercial  quartz  crystal  oscillator  manufacturer's 
approach,  utilizing  low  cost  digital  integrated  circuits  as  basic  oscillator  building  blocks. 
Experience  gained  as  a  result  of  extensive  linear  oscillator  design  work  will  be  applied, 
both  in  the  recommended  packaging  approach,  and  analytical/empirical  design  verifica¬ 
tion. 


2.  DESIGN  APPROACH 

To  combine  performance  with  low  cost,  it  is  necessary  to  synthesize  a  simplistic  oscillator 
design,  utilizing  available  I.C.  technologies.  Previous  experience  with  linear  oscillators 
highlighted  the  need  for  low  power  dissipation  to  aid  integral  mounting  with  the 
instrument.  An  ideal  I.C.  technology  combining  speed  and  low  power  is  CMOS.  The  basic 
oscillator  design  approach  is  to  satisfy  the  Barkhausen  conditions  for  oscillation.  Figure 
12  is  a  block  diagram  illustrating  the  basic  loop  criteria  to  sustain  oscillation  requiring  a 
loop  gain  of  unity  and  zero  degrees  phase  shift.  A  schematic  diagram  of  the  basic 
oscillator  circuit  is  shown  in  Figure  13.  The  crystal  resonator  is  being  operated  in  its 
series  resonant  mode  at  nominally  40  kHz.  A  configuration1  of  the  gated  oscillator  variety 
is  used  whereby  the  open  loop  period  of  oscillation  is  determined  by  the  supply  voltage, 
CMOS  threshold  voltage  and  feedback  RC  time  constant.  The  closed  loop  frequency  is 
determined  by  the  resonator  feedback  element  which,  for  a  high  Q  device,  has  a  large  rate 
of  change  of  phase  W/R  to  frequency.  Initial  calibration  errors  or  environmentaly  induced 
component  drifts  will  therefore  produce  a  small  frequency  pull  to  correct  for  same.  This 
characteristic  will  be  put  to  good  use  in  a  temperature  compensation  circuit  to  be 
described  later. 


3.  THEORY  OF  OPERATION 

The  basic  oscillator  design  is  a  form  of  gated  timing  circuit.  Referring  to  the  schematic 
of  Figure  13,  the  switching  operation  occurs  around  the  threshold  voltage  of  the  logic 
devices  used.  Capacitors  and  C„  form  a  capacitive  divider  network  which  will  apply  a 
portion  of  the  step  voltage  change  at  Vg  to  the  input  of  A...  The  voltage  at  Vf  will 
decrease  with  a  time  constant  proportional  to  R  ..  When  the  logic  threshold  voltage  of 
A. .  is  reached,  both  the  V .  and  V„  will  exhibit  a  reverse  polarity  voltage  step  which  will 
decay  with  the  same  time  constant:  A  diagram  of  this  switching  action  follows. 
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G  (jw)  H(jw)  =  1  +  jO 

FIGURE  12.  BARKHAUSEN'S  CRITERION  FOR  S.S.  OSCILLATION 


R1  R2  R3 


FIGURE  13.  SATURATING  OSCILLATOR  SCHEMATIC 
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The  sum  of  the  switching  periods  T,  +  are  set  with  the  external  trim  resistor  (R  .) 
and  should  be  close  to  the  resonators  series  resonant  period.  The  high  Q  resonator  tuWd 
circuit  will  pull  the  closed  loop  frequency  near  its  series  resonant  point.  The  hybrid  used 
also  contains  a  current  regulator  diode  in  series  with  the  B  line  to  minimize  the  effects 
of  line  voltage  and  circuit  power  consumption  changes  on  the  critical  frequency 
determining  parameters. 

4.  TEMPERATURE  COMPENSATION 

As  described  in  Appendix  E,  the  temperature  sensitivity  of  the  resonator  bias  frequency 
can  be  modeled  as  a  downward  sloping  parabola  centered  around  a  frequency  at  which  the 
temperature  coefficient  is  zero  (i.e.,  Zero  Temperature  Coefficient  (ZTC)  or  turnover 
point).  A  perfect  match  of  two  resonators  with  the  same  ZTC  point  and  the  same  second 
order  coefficients  would  produce  a  difference  frequency  with  a  ZTC  over  the  full  military 
environment.  In  practice  the  second  order  coefficients  of  two  resonators  match  quite 
well  with  compensation  required  for  the  different  ZTC  points.  It  can  be  shown  that 
adding  a  function  which  linearly  varies  with  temperature  to  a  parabolic  function  can  shift 
the  latter's  ZTC  point.  Varying  the  gated  time  constant  resistance  (R  J  linearly  as  a 
function  of  temperature  will  produce  a  linear  change  in  output  oscillator  frequency  with 
temperature  in  the  present  digital  oscillator  configuration  (see  Figure  5-2).  What  remains 
to  be  done  is  to  determine  the  scale  factor  of  this  frequency  vs  temperature 
characteristic  and  to  design  a  resistance  network  to  satisfy  these  conditions.  The  choices 
are  metallic  positive  T/C  resistors  or  a  linearized  version  of  the  negative  temperature 
coefficient  thermistors.  The  latter  configuration  was  chosen  due  to  the  more  compatible 
resistance  values  at  the  temperatures  of  interest. 

In  practice,  an  existing  commercial  hybrid  module  was  used  for  the  basic  oscillator 
section  which,  together  with  accompanying  discrete  trim  components,  was  integrally 
mounted  on  the  low  cost  instruments.  A  schematic  of  this  brassboard  configuration  is 
shown  in  Figure  14.  It  should  be  noted  that  only  one  of  the  two  oscillators  require  the 
linear  temperature  compensation  network.  The  resonator  with  the  lower  ZTC  point  is 
compensated. 

When  driving  this  resonator  with  a  saturating  type  of  oscillator,  in  the  series  resonant 
mode,  the  change  of  frequency  vs  R  .  is  as  follows: 


FREQ 


^ext 


For  a  reasonable  range  of  external  resistance  the  frequency  change  curve  is  linear. 
Without  compensation  the  beam  frequency  may  be  accurately  modeled  as: 
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FIGURE  14.  LOW  COST  VBA  BRASSBOARD  SCHEMATIC 
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f  (T)  =  A0  +  AXT  +  A2T2 

Completing  the  square  of  Equation  (6)  gives: 

-A  2  2 

f  (T)  =  A Q  al  +  A2  (T  +  Al 

4A2  2A2 
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The  inflection  point,  or  ZTC  point  of  Equation  (7)  is  at 


T 

ZTC 


With  the  addition  of  a  linear  resistance  change  vs  temperature,  using  the  R  .  connection, 
a  linear  KT  term  is  added  to  Equation  (6)  giving:  ext 


/  2 

f  (T)  =  Aq  +  AXT  +  A2T  +  KT 


Again  completing  the  square  gives: 

f'  (T)  =  A0  ~  (At  +  K)2 
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Therefore  the  compensated  beam  has  a  ZTC  point  at 


T 

1  ZTC 

A  tZtc 


-  (At  +  K) 
2A2 

-K 
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Thus,  linear  temperature  compensation  has  the  effect  of  shifting  a  given  resonator's  ZTC 
temperature.  It  should  be  noted  that  the  second  order  coefficient,  A„,  is  preserved  by  the 
compensation.  Graphically  the  impact  of  compensation  is  as  follows: 


FREQ 


30 


This  compensation  technique  provides  a  means  for  matching  the  ZTC  temperature  of  two 
separate  resonators.  Only  one  of  the  two  resonators  need  be  compensated  and  providing 
we  have  a  match  of  the  second  order  coefficients  the  difference  frequency  shall  be 
virtually  constant  over  temperature. 

To  provide  this  compensation  a  passive  load  resistance  network  must  be  implemented 
which  varies  as  a  linear  function  of  temperature  and  is  substituted  for  Rext  in  Figure  13. 
A  linearizing  thermistor  configuration  is  as  follows:  Al 


R2 
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The  expression  for  Rgxt  is: 


R 


ext 


ext 


=  +  Rt  //  r2 

Rl  +  Rt  R2 

Rt  +  r2 


(11) 
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Graphically  the  resistance  vs  temperature  characteristic  of  the  network  is: 

^ext 


ROOM  TEMP 


A  large  portion  of  the  curve  centered  about  the  midpoint  (room  temperature)  is  linear. 
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The  thermistor  R^.  is  an  exponential  function  of  temperature  and  can  be  approximated 
with  a  series  expansion  resulting  in  the  following  expression  for  Rgxt: 


R  .  =  R, 
ext  1 


+  R  R-> 

_Q _ 

Ra  +  R- 


+  Rn  R0  •>  a  A  T 


(13) 


where  Ro  =  Room  Temp,  resistance  of  thermistor 

a  =  Temperature  coefficient  of  thermistor 


Let  Ro  f  ~  Ri  +  Ro  r2 

Kq  +  K2 


(14) 


(K0  +  R2)2 


Rewriting  Equation  (13)  gives: 


(15) 


Rext  =«o',kj‘T  «« 

The  above  equation  behaves  as  a  linear  function  of  temperature  over  a  sizeable 
temperature  range. 

5.  EXPERIMENTAL  RESULTS 

The  temperature  compensation  of  accelerometer  bias  using  this  technique  was  demonstra¬ 
ted,  as  described  in  Section  4.  The  procedure  was: 

o  The  temperature  sensitivity  of  the  individual  resonators  of  LC1  and  LC3  was 
determined.  This  data  appears  in  Figure  8  A  and  B  with  the  resultant 
temperture  sensitivity  of  the  difference  frequency  which  appears  in  Figure  8 
C. 

o  Using  this  data  and  the  procedure  described,  the  characteristics  of  the  R 
circuit  was  determined  and  then  installed.  ex 

o  The  temperature  tests  of  Figure  8  were  then  repeated  for  these  new 
conditions.  As  Figures  8  E  and  F  demonstrate,  the  procedure  was  very 
effective  with  the  bias  temperature  sensitivity  being  reduced  to  the  order  of 
1  Pg/°F  over  the  limited  temperature  range  of  the  experiment. 

6.  RECOMMENDED  FUTURE  ACTIVITIES 

Future  activities  should  include  a  further  refinement  of  the  linear  temperature  compensa¬ 
tion  scheme  whereby  the  full  military  temperature  range  (i.e.,  -55°C  to  +71°C)  is 
addressed.  Compliance  over  this  environment  will  make  the  low  cost  VBA  a  prime 
candidate  for  future  strapdown  systems. 


Short  term  stability  testing  of  these  LC1,  LC3  units  revealed  that  package  size  and 
physical  proximity  played  an  important  part  in  the  results.  When  a  discrete  equivalent 
circuit  was  used  for  breadboard  evaluations,  there  was  a  10/1  degradation  in  peak  to  peak 


drift.  One  of  the  major  contributors  to  this  instability  is  humidity  changes  which  will 
affect  interlead  capacitances.  From  analytical  studies  conducted  on  a  linear  oscillator, 
the  resonator  drive  lead  capacitance  had  the  highest  sensitivity  factor.  It  should  be  noted 
that  a  change  in  lead-to-lead  capacitance  has  a  direct  effect  on  the  value  of  the 
resonator's  electrode  capacitance  (i.e.,  Co  =  1.5  pf.  nom).  As  the  value  of  R„  increases 
the  sensitivity  also  increases.  A  solution  to  this  problem  is  to  develop  a  small,  low  cost 
hybrid  module  that  is  mounted  directly  to  the  low  cost  VBA  unit.  With  careful  layout  am- 
short  stiff  interconnect  system  this  hybridized  version  will  optimize  performance.  A 
preliminary  study  on  parabolic  temperature  compensation  is  included  in  Appendix  H.  This 
type  of  compensation  will  essentially  make  each  resonator  insensitive  to  temperature 
variations,  thereby  increasing  yield.  A  cost  and  size  trade-off  study  would  be  conducted 
as  part  of  this  effort. 


SECTION  VI 

VBA  RELIABILITY/MAINTAINABILITY  PREDICTION 

1.  INTRODUCTION 

A  reliability  prediction  was  generated  to  comply  to  para.  4.3  of  the  Statement  Of  Work 
for  the  VBA  Accelerometer  Study  Program.  Maintainability  aspects  of  the  design  are  also 
addressed  herein.  The  Reliability  Prediction  was  performed  for  a  single  beam  acceler¬ 
ometer  for  an  Airborne  Fighter  application.  The  prediction  considers  both  inhabited  and 
uninhabited  flight  environments.  This  prediction  is  based  upon  the  latest  design 
configuration  and  utilizes  the  failure  rates  of  MIL-HDBK-217C-1  and  the  failure  rates 
derived  from  Kearfott  field  performance  data  on  similar  components  in  inertial  instru¬ 
ments.  M1L-HDBK-217C-1  environmental  'K'  factors  were  used  to  adjust  failure  rates  for 
Inhabited  and  Uninhabited  Flight,  where  applicable,  and  MIL-STD-756A  environmental  K- 
factors  for  parts  not  included  in  MIL-HDBK-217C-1. 

Four  different  VBA  ambient  operational  temperatures  were  also  considered  for  this 
analysis. 

2.  RESULTS  AND  CONCLUSIONS 

The  prediction  analysis  results  are  summarized  in  Table  1.  The  failure  rate  for  the 
Uninhabited  Fighter  application  varies  from  5.7047  f/10®  to  6.6357  f/10^  hours  for  the 
four  ambient  temperatures  considered  (55°C,  60°C,  65°C  and  70°C).  This  translates  to  a 
Mean  Time  Between  Failure  (MTBF)  of  175,293  to  150,699  hours  with  increasing 
temperature. 

The  failure  rate  for  Inhabited  Fighter  application  varies  from  5.0667  f/106  to  5.7161  f/106 
hours  with  ambient  temperature  increasing  from  50°C  to  70°C.  This  translates  to  an 
MTBF  of  197,366  to  174,945  hours. 

The  hybrid  oscillator  is  the  predominant  influencing  factor  in  the  change  in  both  the 
failure  rate  and  MTBF  with  temperature  and  environments. 

The  VBA  accelerometer  is  a  sealed  component.  Therefore,  it  has  no  maintainability 
requirements. 

3.  DISCUSSION 

a.  SINGLE  BEAM  ACCELEROMETER 

(1)  Parts 

The  parts  comprising  the  single  beam  accelerometer,  quantity  used,  and  sources  of  the 
failure  rate  data  are  listed  below  in  Table  2.  For  the  purposes  of  this  analysis,  a  military 
equivalent  of  the  hybrid  oscillator  was  considered. 

(2)  Failure  Rates 

Failure  rates  were  obtained  from  MIL-HDBK-217C-1,  and  also  derived  from  Kearfott  field 
performance  data  on  similar  parts,  where  MIL-HDBK-217C-1  failure  rate  data  was  not 
available.  The  failure  rates  of  the  following  VBA  parts  are  based  upon  a  field  evaluation 
of  2,938  Kearfott  Single  Axis  Accelerometers  which  accumulated  6,822,759  equivalent 


laboratory  environment  operating  hours.  These  accelerometers  were  used  on  14  different 
programs  for  missiles,  aircraft,  and  space  applications. 

TABLE  2.  PARTS  AND  (FAILURE  RATE)  SOURCE  OF  DATA 


PART 

QUANTITY 

SOURCE  OF  DATA 
(FAILURE  RATE) 

HYBRID  OSCILLATOR 

1 

MIL-HDBK-2 1 7C- 1 

THERMISTER 

1 

MIL-HDBK-2 1 7C-1 

i  RESISTOR,  COMP.  (R) 

1 

MIL-HDBK-2 1 7C-1 

CRYSTAL 

1 

MIL-HDBK-2 1 7C-1 

CONNECTOR 

1 

MIL-HDBK-2 17C-1 

FLEXURE 

1 

KEARFOTT  DATA 

FEED-THRU 

2 

KEARFOTT  DATA 

SEAL 

4 

KEARFOTT  DATA  | 

] 

PARTS  FAILURE  RATE 

Flexure  0.2565  f/10®  h 

Seal  0.0321  f/lo£  h 

Feed-through  0.0641  f/10b  h 

MIL-STD-756A  environmental  adjustment  'K'  factors  were  applied  to  the  laboratory 
failure  rates  for  these  parts  to  obtain  the  equivalent  fighter  aircraft  environment  failure 
rate. 

(3)  Reliability  Prediction  Calculation 

Reliability  predictions  were  performed  for  airborne  fighter  aircraft  at  four  different  MV  A 
ambient  operational  temperatures  (similar  to  Kearfott  Single  Axis  accelerometer  applica¬ 
tion)  of  55°C,  60°C,  65^,  and  70°C  for  both  Inhabited  Flight  and  Uninhabited  Flight 
environment. 

The  MTBF's  and  failure  rates  for  each  ambient  temperature  have  been  calculated  for  both 
Inhabited  and  Uninhabited  Flight.  The  prediction  details  are  presented  in  Tables  3  and  4, 
respectively. 
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TABLE  3.  ENVIRONMENT  -  AIRBORNE,  INHAB,  FTR 


MTBF  (MEAN  TIME  BETWEEN  FAILURE)  (HOURS)  197,265  190,768  183,287  174,945 


TABLE  4.  ENVIRONMENT  -  AIRBORNE,  UNINHAB,  FTR 


MTBF  {MEAN  TIME  BETWEEN  FAILURE  (HOURS)  175,293  167,869  159,641  150,699 


b.  DUAL  BEAM  ACCELEROMETER 


Based  on  engineering  discussions,  it  has  been  determined  that  the  Dual  Beam  Acceler¬ 
ometer  complexity  will  be  twice  that  of  the  Single  Beam  Accelerometer  which  results  in  a 
failure  rate  twice  that  of  the  single  beam  listed  at  any  particular  ambient  temperature  in 
Tables  3  and  4.  This  results  in  a  MTBF  equal  to  one-half  that  listed  in  the  respective 
tables. 

c.  MAINTAINABILITY 

The  single  and  dual  beam  accelerometers  are  sealed  components;  therefore,  they  have  no 
maintainability  requirements.  The  replaceable  items  would  be  returned  to  the  depot  for 
repair. 
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SECTION  VII 

FUTURE  DEVELOPMENT  AREAS 


In  this  section,  the  recommended  tasks  for  future  development  are  made.  Because  the 
dual  beam  approach  offers  the  best  immediate  potential,  we  recommend  that  future 
efforts  proceed  along  these  lines.  The  individual  areas  of  effort  to  achieve  the  ultimate 
goals  are  described  below. 

1.  GLASS  FRIT 

This  effort  is  intended  to  eliminate  the  use  of  epoxy  in  installing  the  crystal  resonators 
into  the  VBA.  Initial  development  work  has  been  performed  as  part  of  another  program  to 
achieve  a  quartz  to  metal  bond  using  a  low  melting  temperature  glass  (glass  frit).  Sample 
joints  have  been  made  of  quartz  to  monel.  Monel  has  been  found  to  be  a  good  choice  for 
this  joint  because  it  is  a  good  expansion  coefficient  match  to  quartz  and  it  also  forms  the 
proper  oxide  that  promotes  wetting  by  the  melted  glass.  The  monel/quartz  assembly  will 
be  fabricated  in  a  sequence  of  operations  as  illustrated  in  Figure  15.  Figure  15  illustrates 
how  a  low  melting  temperature  glass  will  form  the  joint  between  the  quartz  and  monel. 
The  monel  can  then,  in  turn,  be  welded  to  the  VBA  assembly.  This  joint  will  then  be 
capable  of  withstanding  the  200°C  bakeout  needed  to  achieve  the  bias  stability  goals.  The 
joint  will  not  exhibit  the  plastic-like  behavior  of  the  epoxy,  but  will  result  in  very  stable 
beam  boundary  conditions. 

2.  ALL-WELDED  ASSEMBLY 

Future  VBA  builds  should  be  of  an  all-welded  design  so  that  a  200°C  bakeout  can  be 
achieved  and  organics  can  be  completely  eliminated  from  the  assembly.  It  is  estimated 
that  improved  resonator  mounting,  as  described  in  the  previous  paragraph,  along  with  high 
temperature  bakeout  and  better  resonator  matching,  can  reduce  the  frequency  difference 
trends  to  significantly  less  than  1  fig  per  day.  The  all-welded  design  effort  will  consist  of 
redesigning  the  VBA  assembly  joints  to  eliminate  the  epoxy  seals  and  replace  them  with 
EB  welds. 

3.  GAS  DAMPING 

The  VBA's  evaluated  as  part  of  this  program  are  filled  with  1/20  of  an  atmosphere  of 
neon.  There  is  some  indication  that  helium  or  perhaps  a  mixture  of  helium  and  neon  will 
provide  a  better  combination  of  high  resonator  Q  and  high  proof  mass  damping.  Part  of 
this  effort  will  be  to  further  analyze  the  damping  mechanism  and  perform  various  other 
experiments  to  further  evaluate  the  effects  of  gas  pressure  and  gas  viscosity  on  the 
resonator  and  the  proof  mass  damping  mechanism. 

4.  MINIATURIZATION 

The  envelope  of  the  breadboard  VBA  is  too  large  to  make  the  VBA  interchangeable  with 
currently  available  accelerometers  in  many  applications.  The  present  size  is  dictated  by 
the  present  geometry  of  the  beams  used.  It  is  suggested  that  a  miniaturization  effort  be 
undertaken  with  the  goal  of  achieving  a  dual  beam  VBA  design  including  electronics 
having  an  overall  length  of  about  1  inch  and  a  diameter  of  1  inch  or  less.  This  effort  will 
include  the  design  and  evaluation  of  a  miniature  resonator  followed  by  the  design  of  a 
miniaturized  version  of  the  dual  beam  VBA. 
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5.  HYBRID  OSCILLATOR 


During  the  course  of  development  tests  both  hybrid  and  discrete  component  versions  of 
various  oscillators  were  evaluated.  It  was  discovered  that  the  small  hybrid  oscillators 
performed  at  least  an  order  of  magnitude  better  than  the  discrete  component  versions  of 
the  same  oscillators.  This  performance  improvement  of  hybrids  over  discrete  component 
versions  is  attributed  to  the  reduction  of  stray  capacitance  effects  by  an  order  of 
magnitude  or  more.  For  both  performance  and  envelope  purposes  it  is  suggested  that  a 
small  hybrid  version  of  the  oscillator  design  to  be  selected  be  developed  and  mounted 
within  the  VBA  body. 

6.  RESONATOR  GEOMETRY  IMPROVEMENT 

The  present  resonator  machining  processes  result  in  a  beam  geometry  which  is  not 
completely  symmetrical  in  cross  section  and  also  not  completely  symmetrical  at  the  beam 
root.  It  can  be  shown  that  these  non-symmetrical  conditions  makes  the  resonator  more 
susceptible  to  instabilities  of  the  beam  boundary  conditions.  It  is  therefore  suggested  that 
the  resonator  fabrication  procedure  be  improved  to  result  in  more  favorable  resonator 
geometry. 

7.  FINITE  ELEMENT  ANALYSIS 

Using  finite  element  analysis  has  been  a  great  advantage  in  determining  resonator  static 
and  dynamic  behavior.  Some  preliminary  work  has  also  been  performed  :n  using  finite 
element  analysis  to  also  predict  and  adjust  the  temperature-frequency  characteristics  of 
the  resonator.  Presently  we  are  able  to  analytically  predict  the  shape  of  the  resonator 
frequency  vs.  temperature  characteristics,  but  we  are  not  able  to  predict  the  turnover 
(ZTC)  temperature.  This  is  because  the  stresses  and  strains  at  the  root  of  a  vibrating 
beam  are  very  complicated  and  it  is  in  the  root  where  the  major  portion  of  the  vibration 
strain  energy  is  concentrated.  A  detailed  finite  analysis  program  is  available  which 
considers  the  entire  quartz  compliance  tensor.  It  is  also  capable  of  performing  coordinate 
transformations  of  this  tensor  which  result  when  different  angular  cuts  of  quartz  are  used. 
This  program  is  also  capable  of  running  in  double  precision  which  is  needed  to  resolve 
small  frequency  changes  experienced  due  to  the  temperature  changes  in  the  various  terms 
of  the  compliance  tensor.  It  is  believed  that  by  using  this  analysis  method,  a  prediction  of 
the  turnover  temperature  vs.  crystal  cut  can  be  made.  The  goal  of  this  effort  will  be  to 
increase  the  turnover  temperature  to  a  more  favorable  temperature  and  also  determine  a 
crystal  cut  which  has  the  lowest  temperature-frequency  sensitivity. 

8.  NON-PRISMAL  RESONATOR  GEOMETRY 

To  date  only  prismal  resonator  geometries  have  been  used  since  they  are  the  simplest  to 
analyze  and  to  machine.  However,  the  analysis  capability  and  the  machining  capability  is 
at  hand  to  analyze  and  machine  resonators  with  non  prismal  beams.  There  is  some 
evidence  to  show  that  non  prismal  beams  will  have  a  higher  force-  frequency  sensitivity 
and  also  be  less  sensitive  to  changes  in  resonator  boundary  conditions.  Higher  scale  factor 
and  improved  bias  stability  should  result. 

9.  LINEARITY  IMPROVEMENT 

It  has  been  discovered  that  subtle  changes  in  the  way  the  axial  load  is  applied  to  the 
resonator  can  affect  the  resonator  linearity.  By  design,  a  resonator  loading  mechanism 
can  be  devised  which  will  tend  to  cancel  the  second  order  frequency-force  sensitivity. 


The  net  result  will  be  a  resonator  frequency-force  sensitivity  that  is  much  more  linear 
than  the  present  resonators.  This  improvement  will  make  the  task  of  matching  and 
trimming  individual  resonators  and  proof  masses,  in  order  to  null  out  the  effects  that 
much  easier. 
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LIST  OF  SYMBOLS 


‘'I'  c2 '  c3»  c4 >  Coefficient!  of  beam  equation  determined  from 

boundary  conditions 

M  “  Bending  Moment 

V  ■  Shear  Reaction 

PE  -  Potential  Energy 

KE  -  Kinetic  Energy 

I  ■  Section  Moment  of  Inertia 

p  *  Density 

y  «  Time  derivative  of  y 

m  *  Mass 

Time  derivative  of  8 
ao»  aj_,  a2/  a3  -  Beam  power  series  coefficients 
G  -  Shear  Modulus 

(«.,  t,  b)B  -  Length,  thickness,  and  width  of  beam  (B) 

E  -  Elastic  Modulus  8.61  (lO**)  dyne/cm  for  quartz 

a  m  Beam  vibrational  frequency  in  rad/sec 

N  -  Beam  vibration  mode  (N  -  1  for  fundamental) 

y  ■  Beam  amplitude 

y1  •  dy/dx  (beam  slope) 

y"  -  d2y/dx2  (beam  curvature) 

y"»-  d3y/dx3 

x  -  Location  along  beam  length 

Other  symbols  as  defined  in  various  illustrations  and  text. 


INTRODUCTION 


In  determining  the  vibrational  character istics  of  a  quartz  resonator, 
piezoelectric  effects  are  neglected.  This  is  justified  because 
of  the  high  Q  (about  40,000)  of  the  resonating  system.  The  analysis, 
therefore,  considers  only  Hook's  Law  and  Newton's  Law. 

The  analysis  uses  the  Rayleigh  Method  which  has  been  modified  to 
include  beam  tension  effects.  The  reason  is  that  the 
Rayleigh  Method  is  easy  to  follow  and  it  agrees  well  with  experi¬ 
mental  data  and  more  rigorous  analytical  methods.  It  also  more 
easily  adapts  itself  to  non-prismal  shaoes  and  can  also  be  used 
to  assess  the  effects  of  shear  and  rotary  inertia  which  is  usually 
neglected  in  flexure  beam  analysis. 
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BEAM  ANALYSIS  CONSIDERING  FLEXURE /  TENSION  AND  LINEAR  DISPLACEMENT 


2.1  Procedural 

The  Rayleigh  Method  will  be  used  and  the  basic  steps  are  as 
follows : 

1.  Assume  a  beam  deflection  shape  that  is  reasonable. 

2.  Assume  harmonic  motion. 

3.  Determine  and  then  equate  the  maximum  potential  energy  of 
the  vibrating  system  to  the  maximum  kinetic  energy. 

4.  Solve  for  the  resonant  frequency  of  the  system.  This  can 

be  done  since  the  kinetic  energy  term  will  include  the  resonant 
frequency  term. 

The  analysis  of  this  section  will  consider  only  the  potential  energy 
due  to  flexure  and  tension  and  only  the  kinetic  energy  due  to  linear 
velocity.  The  effects  of  shear  and  rotation  are  neglected  in  the 
analysis  but  are  commented  upon  later. 

2.2  Assumed  Beam  Deflection 

The  assumed  deflection  will  be  that  of  a  fixed-fixed  vibrating  beam 
where  tension  effects  are  considered  to  have  a  negligible  effect 
on  the  deflection  shape.  The  deflection  shape  for  the  first  five 
modes  is  shown  in  Figure  2-1.  The  general  expression  for  these 
shapes  is  given  in  Reference  1  and  also  in  other  vibration  literature 
as  Equation  (2-1). 

CX4uJv.  ( 2_1) 


MODE  NUMBER 


t  1 


(EVEN  MODES) 


FIGURE  2-1 

BEAM  DEFLECTION  CURVES 
FOR  MODES  1  THROUGH  5 
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2.2  (cont'd) 


Prom  odd  and  even  symmetry  of  Figure  2-1  it  can  be  seen  that  the 
coefficients  Ci  and  C3  will  be  zero  of  odd  modes  (even  functions, 
N  ■  1,  3,  5...)  and  the  coefficients  C2  and  C4  will  be  zero  for 

the  even  modes  (odd  functions,  N  -  2,  4,  6...).  The  remaining 

coefficients  are  evaluated  for  the  boundry  conditions  of  zero 
deflection  (y  -  0)  and  zero  slope  (dy/dx  -  0)  at  the  beam  ends 
(x  ■  -1/2) .  Equation  (2-1)  then  reduces  to  Equations  (2-2)  and 


2-3) .  The  values  of  Table  2-1  will  satisfy  the  boundary  conditions 


for  the  Various  modes. 


For  odd  modes  (even  functions) 


(2-2) 


For  even  modes  (odd  functions) 


(2-3) 


TABLE  2-1.  BOUNDARY  CONDITIONS  FOR  MODES  OF  OPERATION 


MODE  MODE  (M) 


FUNDAMENTAL 


1.0037  (3f/24-) 


.1328 


.5999  (  511  /2 lB) 


.0278 


1.0000  (7V2£B) 


-.00580 


1.0000  ( 9 1 / 2  £B ) 


-.00120 


1.0000  (11H/2*b)  -.000250 


I 
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2.2  (cont'd) 

Not#  that  an  approximation  for  k  ia 

iO 

Also  note  that  the  actual  values  of  C3  and  C4  will  depend  on  the 
beam  vibration  amplitude. 

Tables  2-2,  3  and  4  include  the  successive  differentiations  of 
Equations  (2-2  and  3)  with  respect  to  x  (y1,  y"  and  y" ' )  along  with 
values  at  various  locations  along  the  beam  for  modes  1,  2  and  3. 
These  functions  are  plotted  in  Figures  2-2,  3  and  4. 
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TABLE  2-2 
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TABLE  2-3 


MODE  2  CHARACTERISTICS 


FIGURE  2-3.  MODE  2  CHARACTERISTICS 


TABLE  2-4 


MODE  3  CHARACTERISTICS 
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FIGURE  2-4.  MODE  3  CHARACTERISTICS 
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2.3  Flexure  Potential  Energy  (PE)m 


A  email  increment  of  the  deflected  beam  ie  illustrated  in  Figure 


U 


BENDING  MOMENT  EFFECT 
FIGURE  2-5 


The  potential  energy  of  a  system  which  has  a  linear  spring  rate  is 
1/2  the  spring  rate  times  the  square  of  the  deflection.  Therefore, 
the  incremental  potential  energy  of  the  dx  beam  element  is: 


KPeV- 


(2-5) 


From  strength  of  materials  comas  the  following  two  expressions 
whioh  relate  bending  angle  (9)  to  bending  moment  (M)  and  bending 
moment  to  curvature  (yN). 


(2-6,7) 
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2.3  (cont'd) 

Making  the  indicated  substitutions  will  give  the  following 
expression  for  the  total  flexure  potential  energy  for  1/2 


X«0  (2-8)_ 


2.4  Tension  Potential  Energy  (PE)T 

Tension  potential  energy  is  determined  by  considering  the  beam 
as  a  string  under  tension  as  shown  in  Figure  2-6. 


T 


TENSION  EFFECT 
FIGURE  2-6 

The  potential  energy,  in  this  case  the  change  in  potential  energy, 
is  the  tension  force  (T)  times  the  change  in  length  (ds  -  dx)  of 
the  "string"  as  it  deflects  from  its  zero  position  to  its  deflected 
position.  The  incremental  potential  energy  of  the  dx  beam  element 
is : 
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kCKyr-T(4*-A*} 


(2-9) 


2.4  (aont'd) 

From  Figure  2-6  is  derived  the  following  expression  for  ds  along 
with  its  approximation  using  the  binomial  expansion. 


Substituting  Equation  (2-10)  into  (2-9)  results  in  the  following 
expression  for  the  total  tension  potential  energy  for  1/2  the 


beam. 


(pEV  I  J  W* 


(2-11) 


2.5  Linear  Kinetic  Energy  (KE)L 

A  small  increment  of  the  deflected  beam  is  illustrated  in 
Figure  2-7. 


LINEAR  MOMENTUM  EFFECT 
FIGURE  2-7 


60 


2.5  (cont'd) 


The  element  dx  is  at  rest  at  maximum  deflection  and  has  maximum 
velocity  in  the  y  direction  when  y»0.  For  such  a  system  the 
maximum  incremental  kinetic  energy  is  1/2  the  mass  times  the 
square  of  maximum  linear  velocity. 


(2-12) 


For  the  element  dx  experiencing  harmonic  motion: 

UL> 


(2-13,  14) 


2.6 


Making  the  indicated  substitutions  results  in  the  following 
expression  for  the  total  linear  kinetic  energy  for  1/2  the  beam. 


(2-15) 


Determine  Beam  Resonant  Frequency  (u>B) 

As  stated  in  Steps  3  and  4  of  section  2.1,  the  beam  resonant 
frequency  is  determined  by  equating  maximum  kinetic  energy  to 
maximum  potential  energy  and  solving  for  uq  as  follows,  p 
From  E  KE  -  EPE 


(keV* 


Substituting  Equations  (2-C,  11  and  15)  into  (2-16) 
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Solve  for  u 


UJ- 


V'1 

9W  _ 


(2-18) 


2.6  (cont'd) 


Note  that  synunetry  allows  the  limits  of  integration  to  be  over  the 
whole  or  half  the  beam.  For  simplicity,  the  limits  of  integration 
are  omitted. 

The  bias  frequency  (tog)  is  determined  for  the  condition  T  *  0 . 

..In  SuMO* 

•  Lwjssrj 


for  1  -  bt3/12 


taJ? 


Soto. 
5 !*% 


(2-20) 


Equation  (2-18)  now  becomes: 


UJ  =•  UJt, 


/  H  SkjJj-i- 

'  u£  PJC  j-  tfjy. 


(2-21) 
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Using  the  binomial  expansion 


Additional  terms  may  be  obtained  by  further  applying 
the  expansion. 

2.6  (cont'd) 

Equation  (2-22)  can  be  expressed  in  another  form. 

u‘231 

The  aN  coefficients  are  determined  from  Equation  (2-22)  and  the 
integral  values  of  Table  2-5. 

TABLE  2-5 
INTEGRAL  VALUES 


'  / 

'  2 

3 

JY^ 

.2WC?Jt 

'J&TCfjC 

.T-SOCJ-J? 

S 

JJt  cZ/4 

//.(,  cf/jt 

2M  C?/J 

nr  Cp/43 

yercf/jt' 

woct/j? 

6Wc:/r 

/.aWXt/jf 

cn#69<S/jf 

Note:  Limits  of  integrals  are  x-0  to  x*I/2. 


Integrals  were  evaluated  using  Simpson's  rule  on  the 
values  of  Tables  (2-2,  3,  4). 


i-_  Sjtrjv  . . 


To  demonstrate  the  affectivanaas  of  tha  Rayleigh  method^ coefficients ~ 
*0  thru  a3  wara  evaluated  for  Mode  1  and  found  to  be  6.45,  .145, 

.010  and  .0015  respectively .  This  compares  wall  to  the  accepted 
values  of  6.45,  .148,  .012  and  .0018  respectively  that  were 
determined  by  a  very  rigorous  and  hard  to  follow  solution  of  the 
differential  equation  of  a  vibrating  flexure  beam  under  tension. 

This  is  good  agreement  considering  the  graphical  integrations, 
other  coefficients  are  given  in  Table  2-6. 

TABLE  2-6.  COEFFICIENTS  OF  EQUATION  (2-23) 


masjK 

MA  \ 

A, 

/ 

C.+T 

L9(fOA) 

2 

/7.1 

./xx 

-77p) 

/.oQo'1) 

3 

3*1 

.ow 

vx>  - 


2.7 


Typical  Beam  Tension  Coafficianta 

Typical  beam  dimensions  and  quarts  properties  arei 

l  ■  .173  in  ■  .439  cm 

-  .040  in  ■  .102  cm 

t  *  .0052  in  *  .0132  cm 

E  -  8.61  (10*)  dyna/cm2 
P  *2.65  gm/cm^ 


Prom  thaaa  valuaa  the  axpressions  are  evaluated. 


i  '  -  s.*\ox»o^  «u«r' 

'  e 


JL 


For  the  bean  frequency  expressed  as  the  oower  series  of 
F.a aation  (2-2  3)  ,  tvoical  f reauenev-tension  coefficients  are 
given  for  Modes  1#  2  and  3  in  ..able  2—7. 
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TABLE  2-7.  BEAM  FREQUENCY-TENSION  COEFFICIENTS 


COEFF , 

MODE!\. 

NO 

ko 

kHz 

kl 

Hz/dyne 

k2 

Hz/dyne^ 

k3 

Hz/dyne3 

1 

40 

5.7(10"3) 

-4.4  (lO""10) 

6.3  (10~17) 

2 

110 

|  P 

-7. sdo'10) 

9.7 (10-17) 

3 

216 

8 . 5  (10-3) 

-1.5(10“10) 

0 . 6 (10-17) 
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APPENDIX  B 


ENGINEERING 

TECHNICAL 

REPORT 


SQUEEZE  FILM  GAS  DAMPING 
FOR  THE 
LOW  COST  VBA 


BY 

W.  C.  ALBERT 

ABSTRACT 


The  introduction  of  squeeze  film  gas  damping  to 
the  basic  VBA  design  approach  has  the  potential 
of  greatly  simplifying  the  design,  reducing  the 
size,  increasing  the  ruggedness  and  increasing 
the  scale  factor.  This  report  is  an  analysis  of 
the  frequency  sensitivity  effects  of  gas  damping. 
It  is  concluded  that  there  will  be  some  increase 
in  vibration  transmissibility  (up  to  20%)  at 
some  frequencies  but  this  is  considered  tolerable. 
The  analysis  also  shows  that  serious  resonance 
problems  can  be  avoided  as  long  as  the  damping 
corner  frequency  can  be  kept  high  compared  to 
the  undamped  natural  frequency  of  the  system. 

This  frequency  ratio  for  the  design  analyzed 
was  7 . 
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LIST  OF  SYMBOLS 

Xc  =  Case  displacement 

Xjjj  =  Proof  mass  displacement 

m  =  Proof  mass  mass 

KB  =  Resonator  spring  rate 

Kd  =  Gas  film  spring  rate  (frequency  sensitive) 

D0  =  Gas  film  damping  coefficient  neglecting  frequency 
sensitivity 

Dd  =  Gas  film  damping  coefficient  considering  frequency 
sensitivity 

WN  =  Undamped  natural  frequency  neglecting  frequency 
effects 

WNn  _  Undamped  natural  frequency  considering  frequency 
effects 

wc  =  Damping  corner  frequency 

£0  =  Damping  ratio  neglecting  frequency  effects 

Gq  =  Damping  ratio  considering  frequency  effect 

y  =  Gas  viscosity 

P  =  Gas  pressure 

R  =  Damping  gap  radius 

c  =  Damping  gap  (squeeze  film)  thickness 

k  =  wc/wN 

tr  =  Transmissibility  neglecting  frequency  effects 

trQ  =  Transmissibility  considering  frequency  effects 
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1.  INTRODUCTION 

Squeeze  film  gas  damping  is  being  considered  for  the  low  cost  version 
of  the  VBA.  Because  of  gas  compressibility  effects,  the  gas  film  also 
has  spring-like  characteristics.  This  spring  term  increases  with  in¬ 
put  vibration  frequency.  Conversely,  the  damping  characteristics  of 
the  gas  film  decreases  with  increasing  frequency.  This  behavior  of 
the  gas  film  alters  the  response  of  the  VBA  which  would  otherwise  be 
a  simple  second  order  system. 

This  analysis  determines  how  the  transmissibility  of  the  system  to 
case  vibration  is  altered  by  these  effects.  The  following  procedure 
was  used. 

1.  When  the  frequency  sensitivity  effect  on  spring  rate  and  damp¬ 
ing  are  ignored,  the  VBA  becomes  a  simple  passive  second  order 
system.  For  these  conditions  the  damping  gap  physical  charac¬ 
teristics  are  determined  for  a  damping  ratio  of  one  (critical 
damping)  and  a  undamped  natural  frequency  of  300  Hz. 

2.  The  corner  frequency  (3  db  point)  at  which  the  damping  term 
is  turning  into  a  spring  term  is  then  determined. 

3.  Determined  next  is  how  this  corner  frequency  alters  the 
natural  frequency,  damping  coefficient  and  damping  ratio. 

How  these  characteristics  vary  with  input  frequency  and 
corner  frequency  are  plotted  in  Figure  3-2. 

4.  From  the  values  determined  in  Step  3,  the.  transmissibility 
of  the  system  is  determined  and  plotted  in  Figure  3-3. 
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2.  DAMPING  GAP  DESIGN  (Neglecting  Frequency  Effects) 

2.1  Step  1 

The  analysis  of  this  section  will  neglect,  temporarily,  the  fre¬ 
quency  dependence  of  the  damping  and  spring  rate  terms.  The 
spring  rate  and  damping  coefficients  are  constants 


SECOND  ORDER  SYSTEM 
FIGURE  2-1 


and  the  problem  reduces  to  a  simple  second  order  system  as  il¬ 
lustrated  in  Figure  2-1.  For  this  system,  the  undamped  natural 
frequency  and  the  damping  ratio  are  described  by  equations  (2-1) 
and  (2-2)  respectively.  Note  these  are  standard  relationships 
(Reference  2) . 

£1=2.  X'171 


co  = 

N 


N/ 


K 


8 


nm 


(2-1) 


(2-2) 


The  double  squeeze  film  damping  gap  is  shown  schematically  in 
Figure  2-2.  The  expression  for  the  damping  coefficient  of  such 
an  arrangement  (again  temporarily  neglecting  frequency  effects) 
is  described  approximately  by  Equation  2-3  (Reference  1). 


aaMMiteiiiaai 
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R 


sd  —  ->  — 


<SAJ 

£  /  OF  Z) 
VISCOS/T?  yLC, 
PRESSURE  p 


DAMPING  MECHANISM 
FIGURE  2-2 


n  -  7  R. 

u0  ^3 


When  Equation  (2-2  and  3}  are  combined  the  following  expression 


is  obtained. 


3  ^ 


(2-4) 


The  damping  gap  c  is  determined  for  the  following  typical  Low 
Cost  VBA  design  conditions. 

=  300  (10-®)  poise  (Neon) 

R  =  0.9  cm 
m  =  10  gm 

WN=  300  Hz  (1.88  (103)  rad/sec) 


co-  1 


The  result  is  a  c  dimension  of  3.61  (10“3)cm  =  1.42  (10~3)  inches 


This  is  considered  a  reasonable  gap  for  manufacturing  purposes. 
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2.2  Step  2 

The  expression  for  the  corner  frequency  in  terms  of  gas  proper¬ 
ties  and  gap  dimensions  is  given  by  Equation  (30)  (Reference  1) 


60  = 


(2-5) 


°  }.07 /ofC 

For  the  typical  conditions  of  Section  2  and  a  pressure  of  1/2 
atmosphere  (P  =  5(10^)  dyne/cm2),  Wc  =  1.3(104)  rad/sec  =  2(103)  Hz. 
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frequency  effects  of  squeeze  film  damping 
3.1  Step  3 

It  is  characteristic  of  a  squeeze  film  gap  to  behave  as  a  pure 
damper  at  low  frequencies  and  then  transform  into  a  pure  spring 
at  high  frequencies.  This  section  will  include  these  effects 
by  modifying  the  simple  system  of  Figure  2-1  into  that  of 
Figure  3-1 


i 


f \Z\j^77" ^ 


'k 


i- 


////'//  A 


Ofr 


Case 

Vibration 


GAS  SPRING  EFFECTS 
FIGURE  3-1 


] 
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Dd  is  the  damping  coefficient  term  which  also  includes  the  fre¬ 
quency  effects.  It  is  described  by  Equation  3-1  (Reference  1) 

Do _ 


"Gff 


(3-1) 


Kq  is  the  spring  effect  of  the  gas  compressibility  which  is  de¬ 
scribed  by  Equation  3-2  (Reference  1) . 

n  ^ 

_  Uo  04 


At  a  particular  input  vibration  frequency,  the  system  of  Figure 


(3-2) 


3-1  will  behave  as  if  it  had  an  overall  spring  rate  and  damping 
coefficient  determined  at  that  particular  frequency.  This  apparent 
natural  frequency  is : 


<^C 


- 

vi  ^ 


(3-3) 


and  after  Equation  (2-1,  2)  and  (3-2)  are  substituted  into  (3-3) 
and  rearranged: 


04 


MO 


(3-4) 


Letting  Wc  -  kw^,  Equation  (3-1)  and  (3-4)  become  respectively: 


(3-5) 


—  *  -UJi if 


UJh) 


o  __ 


/  + 


2  'foJz 


(3-6) 
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Both  these  expressions  are  plotted  in  Figure  3-2, 


The  apparent  damping  ratio  will  also  change  with  input  frequency. 


The  expression  used  for  damping  ratio  are: 


3L'>»com. 


Q-Q 


2  m  wNi 


(3-7) 


Therefore : 


(3-8) 


Equation  (3-8)  is  also  plotted  in  Figure  3-2, 


The  effects  of  DD,  KQ  and  Wc  on  the  transmissibility  of  the  system 
of  Figure  3-1  will  be  determined  next. 

3.2  Step  4 

From  vibration  analysis  (Reference  2)  the  expression  for  the 
transmissibility  of  the  system  of  Figure  2-1  is 

I  -Llktf’  kT  _  2k  (3-9) 

0  I-  .  ,TA  ,  .L  *  -V 


This  expression  is  next  modified  to  account  for  the  frequency 
sensitivity  of  the  damping  and  spring  coefficients  of  Figure  3-1. 


■c.-l  /jO 

'Ui 


(3-10) 


This  expression  was  evaluated  and  plotted  for  =  1  and  various 
values  of  k  in  Figure  3-3. 
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4.  CONCLUS IONS 

From  the  equations  developed  and  the  plots  of  Figures  3-2  and  3-3, 
it  can  be  seen  that  the  frequency  sensitivity  of  squeeze  film  damp¬ 
ing  has  the  following  effect: 

1.  The  apparent  undamped  natural  frequency  increases  with 
increasing  input  vibrational  frequency.  This  is  due  to 
the  squeeze  film  acting  as  a  spring  at  high  frequencies. 

2.  The  damping  coefficient  decreases  with  increasing  fre¬ 
quency. 

3.  The  apparent  damping  ratio  decreases  with  increasing  fre¬ 
quency  . 

4.  If  Wc  is  near  WN,  The  system  will  behave  as  an  underdamped 
system  even  though  it  is  designed  for  critical  damping  at 
low  frequencies.  If  Wc  is  large  compared  to  WN,  then  the 
system  behaves  more  like  a  critically  damped  system  although 
the  resonant  frequency  will  be  higher.  For  the  proposed 
squeeze  film  damped  VBA,  the  ratio  of  Wc  to  WN  is  about  7 
and  therefore  the  transmissibility  will  never  exceed  1.2 

as  illustrated  by  Figure  3~3.  Figure  3-3  also  shows  the 
responce  of  the  system  neglecting  the  frequency  sensitive 
effects. 
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APPENDIX  D 

DUAL  BEAM  VBA  FOUR  POSITION  TEST 


Figure  D-l  is  a  schematic  representation  of  the  Dual  Beam  VBA  test  arrangement.  For  a 
four  position  test  just  bias,  scale  factor,  nonlinearity  and  misalignment  are  considered. 
Therefore,  for  Beam  1,  the  assumed  output  model  and  input  acceleration  are, 
respectively: 

fl  =  f01  +  kllS  +  k2li?2  (D-!> 

and 

g  =  G  sin  (0+0M1)  (D.2) 

Similarly,  for  Beam  2 

f2  =  f02  +  k12  8  +  k22  S2  <D’3> 

g  =  Gsin(0+ 18O°+0M2)  (D.4) 

Where 


f01  and  f02 
kll  andki2 


k2i  and  k22 
0M1  and0M2 


Bias  frequencies  of  Beams  1  and  2,  respectively,  in  Hz 

First  order  scale  factors  of  Beams  1  and  2,  respectively,  in 
Hz/g 

9 

Nonlinearity  of  beams  1  and  2,  respectively,  in  Hz/g 

Misalignment  of  VBA  halves  1  and  2,  respectively,  in 
radians. 


Readings  are  taken  at  6=  0°,  90°,  180°  and  270°.  It  is  assumed  that  misalignment  effects 
on  the  effects  are  a  secondary  effect  on  a  secondary  effect  and  therefore  negligible. 
It  is  also  assumed  that  the  misalignment  effects  are  negligible  for  the  6 =  90°  and  Q-  270° 
positions.  The  model  now  reduces  to  the  following  for  Beam  1  at  the  four  positions: 


<fl>0  =  f01  +  kll  G  6  Ml 

(D.5) 

^^90  =  f01  +  kll  G  +  k21  G 

(D.6) 

(Vl80  =  f01  “k  11  G  6  Ml 

(D.7) 

<fl*270  =  f01  ”  kll  G  +  k21  G 

(D.8) 

A  similar  set  of  equations  is  also  obtained  for  Beam  2. 


By  summing  and  differencing  Equations  (D.5)  through  (D.8),  the  following  expression  for 
the  various  coefficients  are  obtained: 


(D.10) 


vll 


k2 


Ml 


(f.)  '(f,) 

90  270 


(fl> 


90 


Tf— 

+  (fj) 


-  2  f 


270 


01 


2  g 

(fl)  -(fj) 

0  180 


(fi>  -V 

90  270 


(D.ll) 


(D.12) 


The  Dual  Beam  VBA  output  is  the  difference  frequency  of  the  two  beams,  as  described  in 
Section  2.  For  the  difference  frequency  the  total  VBA  output  model  becomes: 


f  =  K0»K,g*K2g 


Where 


Bias  =  Kq  ~  fQi 


f02  in  HZ 


Scale  Factor  =  =  k^  +  k ^  in  Hz/g 

2 

Nonlinearity  =  =  k21  -  k22  in  Hz/g 

The  input  now  is  expressed 
0M1  + 


g  =  G  sin  \d  + 


in  \e* 

2 


(D.13) 

(D.14) 

(D.15) 

(D.16) 

(D.17) 


by  Equation  (D.17),  which  indicates  that  the  misalignment  effects  become  the  average  of 
the  individual  0^  term. 

From  Equations  (D.14)  through  (D.17),  the  following  becomes  apparent: 

•  The  bias  is  nominally  zero  and  the  tolerance  on  the  individual  f^  mismatches 
will  determine  its  actual  value.  The  important  feature  of  common  mode 
rejection  is  also  achieved. 

•  The  scale  factor  is  nominally  double  that  of  the  individual  resonator. 

•  The  nonlinearity  is  also  nominally  zero  and  the  tolerance  on  the  individual  k2 
terms  will  determine  the  actual  value.  Additional  zero  trim  of  this  term  can 
also  be  achieved  by  changing  slightly  the  mass  of  one  of  the  proof  masses. 

•  The  net  misalignment  is  the  average  of  the  individual  misalignments.  These 
terms  can  also  be  adjusted  to  zero  by  conventional  means. 
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APPENDIX  E 

BIAS  TEMPERATURE  SENSITIVITY  MODEL 

For  the  dual  beam  VBA,  a  low  bias  temperature  coefficient  is  achieved  by  matching  the 
individual  beam  temperature  sensitivities.  The  temperature  coefficient  of  the  proposed 
resonator,  like  that  of  all  crystal  resonators,  is  nonlinear.  The  proposed  resonators  have 
an  empirically  determined  frequency  vs  temperature  characteristic  as  illustrated  in 
Figure  E-l. 


f  (FREQUENCY) 


RESONATOR  BIAS  FREQUENCY  VS  TEMPERATURE 
FIGURE  E-l 

The  curve  is  described  very  well  by  a  second  order  expression  having  the  following  typical 
form: 


f  =  f  +  kT  T2  (E.l) 

o  r  z 

The  following  definitions  now  apply  to  Figure  E-l  and  Equation  (E.l). 

TZtC  is  the  Zero  Temperature  Coefficient  (ZTC)  temperature,  sometimes  called  the 
turnover  temperature,  at  which  the  resonator  becomes  virtually  temperature  insensitive. 
The  temperature  at  which  turnover  occurs  is  adjustable  to  some  extent  by  the  proper 
selection  of  beam  geometry  and  crystal  axes  rotation.  Controlling  the  location  of  this 
point  is  a  major  area  of  investigation  on  the  High  Accuracy  VBA  Program.  It  is  a 
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temperature-controlled  device  and  is  being  designed  so  that  its  operating  temperature 
coincides  with  the  ZTC  temperature.  This  will  make  the  High  Accuracy  VBA  virtually 
temperature  insensitive  for  small  temperature  variations. 

AT™  is  the  mismatch  of  individual  resonator  ZTC  temperatures  for  the  dual  beam  VBA. 
As  tne  following  analysis  will  show,  it  is  convenient  to  use  and  T„2  as  temperature 
differences  to  the  accelerometer  operating  temperature  T^. 

The  VBA  bias  is  the  difference  frequency  f^  -  f2  -  Af  and  so  from  Figure  E-l  considering 
just  bias  and  bias  thermal  effects 

Af  =  f  -  f  =  f  -f  +k  -  u  T2 

*1  12  *01  l02  KT1  1  Z1  kT2  1  Z2  (E.2) 


T„j  +  AT ™  is  substituted  for  T Z2  and  after  expansion  and  rearranging,  Equation  (E.3)  is 
obtained. 

Af  =  f01  ~  f02  +  (kTl“kT2^  T  Z1  ~  2kT2  TZ1  ATZM  +  kT2  ATZM  *E'3) 


A  change  in  accelerometer  temperature  T^  will  esentially  be  a  change  in  T^  since 
(TzTC^l  *s  a  constant  an^  80  Equation  (E.3)  is  differentiated  with  respect  to  T™ 


_  d(Af)  2(j<  _  k  _  2k  AT 

akT1  T2nZl  ZKT2a  ZM 

dTZl  dlA 


(E.4) 


To  obtain  the  thermal  sensitivity  in  terms  of  equivalent  input  acceleration  (a)  Equation 
(E.4)  is  divided  by  the  accelerometer  scale  factor  (k^  +  k22  -2kj)  (see  Equation  2.3). 

=  — -T.l  I  kT2)TZl  ~  k12^TZM  (E,5) 

Of  the  two  remaining  terms  of  Equation  (E.5),  the  first  is  the  result  of  a  mismatch 
between  resonators  and  the  second  is  the  result  of  the  ZTC  mismatch.  Because  quartz  is 
a  well-ordered  crystal  material,  the  thermal  characteristics  are  very  consistent.  An 
estimate  of  expected  thermal  sensitivity  will  be  made  for  the  following  typical  conditions: 

k,p  =  -1.3(10  3)  Hz/°C2.  This  is  an  experimentally  determined  coefficient. 


kj  =  100  Hz/g.  This  is  typical  of  the  proposed  scaling. 

(k,^  -  krpj  =  0.01  k,p.  It  is  estimated  that  this  coefficient  can  be  matched  to  1%  by 
proper  selection. 


AT™  =  0.25°C.  It  is  estimated  that  a  ZTC  match  of  0.25°C  can  be  made  by  proper 
selection. 


For  a  temperature  controlled  accelerometer  application,  the  ZTC  point  and  the  operating 
temperature  are  designed  to  be  coincident  so  that  T„,  =  0.  For  this  condition  only  the 
AT,™  term  of  Equation  (E.5)  remains  and  for  the  above  typical  conditions,  Aa/AT  =  - 
3.2fig/°C. 

For  no  thermal  control,  the  ZTC  point  is  selected  to  be  at  the  midpoint  of  the  operating 
temperature  range,  so  that  the  expected  thermal  deviation  from  this  point  is  typically 
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±30°C.  For  the  above  typical  conditions  the  T„.  (thermal  error)  term  of  Equation  (E.5) 
varies  from  0  to  +3.9  figrC  at  the  extremes  oiihe  operating  temperature  range  (0  and 
60°C).  This  thermal  error  source,  coupled  with  the  ATZM  term,  indicates  that  the 
expected  thermal  sensitivity  on  bias  will  vary  typically  fromlsero  to  -3.2  +3.9  =  +  0.6  to  - 
7.1  ng/°C  over  the  operating  temperature  range  of  0  to  60°C. 

Because  of  the  characteristics  of  quartz,  the  thermal  error  on  bias  will  be  very 
predictable  and  very  repeatable  for  effective  system  modeling.  To  achieve  the  ultimate 
goal  of  ljig/C,  a  combination  of  resonator  matching  and  system  modeling  is  proposed. 
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APPENDIX  F 
DUAL  BEAM  VBA 
CROSS  AXIS  SENSITIVITY 


Cross-axis  sensitivity  is  caused  by  a  deflection  of  the  proof  mass  due  to  mg.  forces  as 
illustrated  by  Figure  F-l.  The  T  forces  shown  in  the  free  body  diagrams  are  Ihe  tension 
forces  of  the  resonators  represented  by  springs  K.  The  F^  forces  are  those  of  the  flexure. 

From  summation  of  forces  in  the  input  (i)  direction  the  following  expression  for  Tj  and  T2 
are  obtained. 


T1  =  M1  *i  -  FHi 

(F.l) 

T2  =  _M1  gi  '  FHi 

(F  .2) 

From  the  free  body  diagrams  of  Figure  F-l,  a  summation  of  forces  in  the  x  direction, 
along  with  a  summation  of  moments  will  indicate  that 


FHi  =  Mgx  e  -  MgK 


=  U2  gx  gj 

jt  k  ~£\T~ 


(F.3) 


Considering  just  the  Tj  and  T„  effects  on  the  beam  frequency,  and  also  considering  that 
the  dual  beam  VBA  output  is  the  difference  of  the  two  resonator  frequencies,  the 
following  expression  for  the  difference  frequency  is  obtained. 


Af  =  fl -f2  =  kll  Vk12T2 


'  F-4) 


Where  kj  and  k^2  are  the  first  order  frequency-tension  sensitivities  of  resonators  1  and  2 
respectively. 

After  combining  Equations  (F.l)  through  (F.4)  the  following  is  obtained 


Af  =  k 


11 


/  Ml2?igx\  „  /~M 2  gi  gx  -  M22  gj  gx\ 

^‘gr  k/  )  12\  k.«  k /  j 


For  the  k^*  k12  fck^  and  Mg'S:  M ^  Equation  (F.5)  reduces  to 

•2  JVI2.  , 

)  -  u-r-i  I  cr.  l 


4f  =  2k1MBi-k1  S|K, 


(F.5) 


(F.6) 


Note  that  the  gj  gx  effects  tend  to  cancel  (common  mode  reject). 

To  obtain  the  difference  frequency  in  terms  of  indicated  acceleration,  Equation  (F.6)  is 
divided  by  2  k^M  to  yield  the  following 


’ind 


TO  “  '  1/2 


,M 


*i  g„ 


Or  considering  just  the  g.  gx  effects 


*7 

For  the  topical  values  of  M  =  9  g,  K  =  8  (10  )  dyne /cm  and  j£  =  1  cm, 
0.06  Mg/gZ  per  1%  mismatch  of  the  (M/f {£)  terms. 

The  (M/k£)  mismatch  will  be  10%  or  less  and  so  no  more  than  a  0.6  /xg/g2  g 
expected. 


(F.7) 

(F.8) 

^ind^i^x 
xgj  error  is 
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APPENDIX  G 

LOW  COST  VBA  -  PARABOLIC  TEMPERATURE  COMPENSATION 

The  frequency  versus  temperature  characteristic  of  a  typical  beam  possesses  a  parabolic 
signature 


FREQ. 


TEMP 


For  a  typical  beam  (Reference:  "Saturating  Oscillator  Tests",  an  Engineering  Technical 
Report  by  W.  Albert),  frequency  may  be  related  to  supply  voltage  (provided  the  Constant 
Current  Diode,  CCD,  is  shorted) 


Therefore,  an  effective  temperature  compensation  for  a  single  beam  may  be  realized  by 
parabolically  varying  the  supply  voltage  with  temperature 


TEMP 
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After  compensation  the  frequency  versus  temperature  characteristic  becomes 

FREQ. 


-J - ! _ TEMP. 

TZTC 

To  compensate,  electronics  must  be  designed  which  generate  a  beam  supply  voltage  that 
varies  as  a  parabolic  function  of  temperature.  In  addition,  the  Turnover  Temperature 
(Tztc)  and  Second  Order  Coefficient  (SOC)  of  the  electronics  must  be  adjustable  in  order 
to  cancel  a  particular  beam's  frequency  versus  temperature  characteristic. 

Consider  the  product  of  two  linear  functions  of  temperature: 


Let  G2  (T)  —  Gq2  +■  ^  T 

and  G2  (T)  =  Gq2  +  KG2 


So,  the  product  is 

G^  (T)  G2  (T)  =  (Goi  +  KG1 A  T)  (G02  +  KG2^T) 

=  Gqi  Gq2  +  (G0lKG2  +  G02KGi)AT  +  KG1  Kg2 
=  Ao  +  A*A  T  +  A2  (AT)2 

Where  Ao  =  Gqi  Gq2 

Ai  =  GOI  Kq2  +  G02  Kgi 

a2  =  Kci  KG2 


Note  that  a  parabolic  temperature  characteristic  is  obtained.  If  one  completes  the  square 
of  the  right-hand  side  of  the  above  equation,  the  following  expression  results: 


By  inspection, 


Turnover  Temp.  *  T2TC  “ 

Second  Order  Coeff.  «  SOC  *  A2 


The  following  circuit  is  based  on  cascading  linearly  temperature  sensitive  voltage  dividers 
created  from  conventional  NTC  thermistors  to  provide  a  beam  supply  voltage  that  varies 
parabolically  with  respect  to  temperature. 


tztc 

adjust 


VR  (REFERENCE  VOLTAGE) 


Where  - 


G,  (T)  -  G0  +  KG^T^ 

l  Thermistor  Divider  Gains 

g2  (T)  -  gc  -  kgat J 
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Determine  the  beam  supply  voltage  V8  (T)  - 
Gi  (T)  =  Go  +  KGAT 
VX  (T)  =  Gx  (T)  Vr 

=  Go  Vr  =  Kq  vrat 

V2  (T)  =  KL  Vr 

V3  (T)  =  >5  7!  (T)  +^V2  (T) 

"  %  G0  VR  +  %  Kg  VrAT  +  3j  KX  Vr 

■  \  G0  Vr  +  35  Kx  VR  +  Jj  KG  VR  A  T 

G2  (T)  =  Gc  -  KGAT 

V4  (T)  =  G2  (T)  V3  (T) 

=  (G0  -  KgAT)  (*sGq  Vr  +  JjKjVr  +  Js  KgVr4T) 

=  >5Vr  (Gq  -  KqAT)  (Gq  +  K3  +  Kq  A  T) 

"  *VR  <Go  <Go  +  Kl>  +  (goKG  “  GoKG  -  KjKqMT  -  KG2  (AT)2) 

=  ’sVr  (Gq  (G0  +  Kl)  -  KiKGAT  -  Kq2  (AT) 2  ) 

V5  (T)  =  K3Vr 

Vs  (T)  =  -K2  V4  (T)  +  (K2  +  1)  V5  (T) 

=  -*5  K2  Vr^Gq  (Gq  +  KX)  -  KjKGAT  -  KG2  (AT)  +  (K2+1)K3VR 
=  [K3  (K2  +  1)  -%  K2  Gq  (Gq  +  KiT]  Vr 
+  (35  Kl  K2  Kg  Vr)  A  T 
+  {35  K2  Kq2  Vr)  (AD2 

Vg  (T)  =  Aq  +  A^A  T  +  A2  (AT)2 

Where  Aq  =  (k3  (K2  +  1)  -  h  K2Gd  (Gq  +  KX)J  VR 
Al  =  35Kl  K2  Kq  Vr 

A2  =  >5K2  Kg2  Vr 

Note:  a  T  represents  temperature  deviation  about  Room  Temp.  (+25°C) . 


G01  represents  divider  Gi ' s  output  at  Room  Temp. 


The  above  circuit  features  simple  adjustment  of  the  Turnover  Temperature,  Second  Order 
Coefficient  and  supply  voltage  at  room  temperature. 

For  example 

Turnover  Temp.  *  TztC  “A1  *  -At  »  -  K?KgVp 

“SaJ”  2  t%K2KG2vR) 


Assuming  Kq  is  fixed,  the  Turnover  Temperature  may  be  adjusted  simply  by  varying 
divider  Kj. 

Also, 

Second  Order  Coeff.  =  SOC  -  A2  -  %K2Kq^Vr 


Assuming  Kq  and  Vr  are  fixed,  the  Second  Order  Coefficient  may  be  adjusted  by  varying 
resistor  K2R. 

Finally, 

Supply  Voltage  (@Room  Temp.)  =  Vn  (at  =  0  )  =  A 

s  o 

=  (K3  (K2  +  1)  -%K2  Gd  (G0  +  Klj]  VR 


Assuming  G0  and  Vr  are  fixed  and  and  K2  have  been  determined,  then  the  supply 
voltage  may  be  adjusted  by  varying  divider  K3. 

Design  Example 

Assume  the  following  fixed  parameter  values  in  the  circuit 

G0  =  0.5  (Divider  Output  at  Room  Temp  (+25°C) 

Kg  =  0.010°C"1  (Corresponds  to  -4%/°C  Negative  Temp.  Coeff. 

Thermistor  Divider) 

=  0.0056°F"1 
Vr  =  10.0  VDC 


98 


Assume  the  following  parameters  for  the  beam  to  be  compensated: 


tZTC  “AT  =  -25°C  (Where  C,  T  represents  Temp.  Deviation  aoout 

Room  Temp.  (+25°C) 


SOC  =  -12  PPM/°F2 
fQ  =  40,000  Hz 

&f  =  +2.9  Hz_ 

AV 

VDC 

Vs  =  5.0  VDC 


(Supply  sensitivity  obtained  from  "Saturating 
Oscillator  Tests"  Report) 

(Supply  voltage  desired  at  Room  Temp.  (+25°C) 


Since,  TZTC  =  -  At  =  -  K-i  ,  thus, 

2A2  3k g 

ki  =  ~2tztc  kg 

-  -2  (-25°C)  (0 . 010°C“1) 

Kx  =  0.500  V/V 

Matching  the  Second  Order  Coefficients  for  the  beam  and  supply  voltage  temperature 
characteristics  via  the  supply  sensitivity  At/Av  yields, 


(SOC)  (fO) 


A  f  4v 


-A2  =  ~*5K2  kg2  Vr 


Thus,  K2  =  ~  2  (SOC)  (fo) 

(Af/AV)  (K  G2)  (VR) 

=  -  2 (~12PPB/°F2)  (40,000Hz). 

(+2.9  Hz/VDC)  (0.0056UF-1)2  (10.0VDC) 


K2  =  1.056  -h-/yu 

Finally,  Vg  -  [k3  (K2  +  1)  -H  K2  Gq  (G0  +  m)]  VR 
Vs 

Vr  +  H  K2  Gq  (Gq  +  K].)  =  K3  (K2  +  1) 


Thus,  K3 


1 

k2  +  1 
1 

1.056  +  1 


(f^-  +  K2 

C  5 . 0  VDC  +  h 

10 .  o'  vdc 


G0  (Gc  +  Kl)J 
(1.056)  (0.5) 


1 


(0.5  4  0.500) 


K3  *  0.372  V/V 
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